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PREFACE 


Theoretical and experimental results of high-efficiency oscillations 
in an avalanche diode are reported. Detailed theoretical analyses of 
three different modes or types of high-efficiency oscillation in a PIN 
diode are presented. The results for TRAPATT oscillations in a PIN 
diode are the most extensive. 

The characteristics of TRAPATT oscillations in a PIN diode are 
discussed and an approximate semianalytical solution for the diode vol- 
tage waveform is derived when the diode current is a squarewave. It is 
shown that a traveling avalanche zone is not necessary to generate a 
dense "trapped” plasma and that the^ boundary conditions prevent the 
trapped plasma from completely filling the depletion layer. Typical 
voltage waveforms and corresponding diode power , efficiencies, and 
impedances at the fundamental and higher harmonics are presented. When 
the diode current is a square wave, the diode does not necessarily 
exhibit a negative resistance at all higher harmonics. 

A computer program for TRAPATT oscillations in a PIN diode is 
described. Its running time is two or three orders of magnitude less 
than more exact time-domain computer analyses. Typical results of diode 
power, dc-to-RF conversion efficiency, and required circuit impedances 
are presented for several different current waveforms, which are com- 
posed of up to the seventh harmonic of a square wave and the first two 
harmonics of a half-wave sine wave. It is shown that high-efficiency 
oscillations are possible with diode currents composed of only the 
fundamental or the fundamental and one higher harmonic. 
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A second high-efficiency mode is discussed. It is distinguished 
from the ordinary TRAPATT mode by the fact that, at the beginning of a 
cycle, the diode voltage is above avalanche and the diode current is 
small or assumed to be zero, whereas in the ordinary TRAPATT mode the 
diode voltage is initially high but below avalanche and the diode 
current is large. A semianalytical solution for this second high- 
efficiency mode in a PIN diode is derived, assuming a rectangular 
current waveform. No calculations of diode operating parameters are 
reported for this mode. This mode is believed to be more representative 
of the high-efficiency oscillations observed experimentally than is the 
ordinary TRAPATT mode, when the fundamental frequency is 1/2, 1/3, or 
1/4 of the diode transit time frequency. 

The third mode of high-efficiency oscillations discussed in this 
report is similar to the ordinary TRAPATT mode, except that avalanching 
occurs throughout a major portion of a half cycle rather than only once 
at the beginning of a cycle to form the "trapped plasma." A quasi- 
static approximation is employed to derive a semianalytical solution 
for the voltage across the diode during the multiple avalanche portion 
of the cycle. This result is then pieced together with the solution 
for the other portions of a cycle that are similar to the ordinary 
TRAPATT mode to obtain analytical expressions for the voltage waveform. 
The current waveform is assumed to be essentially a square wave except 
the transitions from maximum current to minimum current are not 
instantaneous . 

Preliminary calculations for this mode in a PIN diode indicate 
that the power increases proportionally to the magnitude of the drive 
current with a small decrease in efficiency. Since this mode is 
characterized by large external currents, the optimum circuit impedance 
for the diode will be less than for the ordinary TRAPATT mode. 
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An analysis is also given for a PIN diode, where it is assumed 
that the ionization coefficient is a constant above a critical E-field 
and is zero for values of E less than the critical field. This approxi- 
mate analysis permits one to draw possible patterns of avalanche regions 
in the depletion layer as a function of time. A rule governing admissible 
patterns is derived; use of this rule has shown that many suggested 
patterns are inadmissible. One useful family of admissible patterns 
has been found. It is further shown that, for a given admissible 
pattern, there are some portions of a cycle where there is a great deal 
of freedom in the waveshape of the external current. In these time 
regions the current can be chosen to optimize the dc-to-RF conversion 
efficiency. 

Preliminary experimental results on the high-efficiency oscillations 
are presented and discussed. Two different experimental circuits, which 
used channel -dropping filters to provide independent harmonic tuning, 
are described. These circuits proved to be inadequate for high- 
efficiency operation because of excessive circuit losses. Simpler cir- 
cuits used to produce high-efficiency oscillations are discussed. The 
effects of harmonic tuning in these circuits were determined indirectly 
from impedance measurements. Inexpensive Fairchild FD300 diodes were 
used and efficiencies comparable with those obtained by Chaffin and 
EerNisse were observed. The impedance measurements made on these cir- 
cuits indicated that these diodes wanted to see a circuit impedance 
with a small real part at all higher harmonics. These results do not 
correspond with the theoretical cases considered in the report, where 
the impedance presented to the diode should be infinite at some of the 
higher harmonics. 

These results tend to substantiate that there are several possible 
voltage and current waveforms which all give relatively high efficiencies. 
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Some circuits required for high-efficiency operation are more easily 
realized than are others. 

The authors wish to thank A. I. Grayzel for many stimulating dis- 
cussions and helpful suggestions during the initial phases of this 
work. 
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1. INTRODUCTION 


The purpose of this contract was to determine theoretically the 

operating characteristics of an avalanche diode oscillator when several 

currents, harmonically related, are flowing in the diode. Emphasis was 

1-4 

placed on the high-efficiency modes of oscillation. A secondary 

purpose was to design, construct, and test an avalanche diode experi- 
mental circuit that would permit independent tuning of at least two 
harmonically related frequencies. 

Most of the objectives of this contract were accomplished. Signi- 

3,4 

ficaut theoretical results were obtained for the TRAPATT mode of 

oscillation in a PIN diode. A PIN diode was considered because it was 

felt that the high-efficiency modes were not critically dependent upon 

the depletion layer doping distribution and the PIN diode simplified the 

analyses. Indeed it is shown that TRAPATT oscillation can exist in PIN 

diodes with voltage and current waveforms similar to those calculated 
3,4 + + 

by others for P NN structures. Two other modes or types of high- 

efficiency oscillations are discussed in this report. It appears that 

these modes have not been discussed previously in the literature. 

Approximate semianalyt ical solutions are given for each of these modes. 

Whereas considerable calculations were made for various operating 

characteristics of the TRAPATT mode only preliminary calculations were 

made for these other high-efficiency modes. Two technical papers based 

on the theoretical results of this report have already been accepted 

16 IV 

for publication and a third one is in preparation. 9 

The results of the experimental circuits originally intended in 
this contract were inconclusive. It was hoped that by proper filtering, 
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a few of the low-order harmonics could be tuned independently to optimize 
the efficiency. The circuits constructed met the design goal of pro- 
viding independent tuning but circuit losses and an improper bias supply 
prevented high-efficiency oscillations from occurring. Simpler circuits 
were also constructed in which several harmonic currents could flow but 
could not be controlled independently. The effects of the higher har- 
monics were determined indirectly by measuring the circuit impedances 
as viewed from the diode terminals. These experimental results indicate 
that the observed oscillations do not correspond to the same sets of 
voltage and current waveforms considered theoretically. The difference 
between theory and experiment tends to substantiate the fact that there 
are several distinct operating conditions (hence, different circuit 
configurations) that lead to high-efficiency operation. Also some 
circuit impedances are easier to realize than are others when multiple 
harmonics are involved. Specifically those circuits with a zero real 
part at higher harmonics are easier to build at microwave frequencies 
than are those circuits that are open-circuited at higher harmonics. 

This report is organized into eight sections including the intro- 
duction, conclusion, and recommendations. In the second section we 
derive a semianalytical solution for TRAPATT oscillations in a PIN diode. 
Equal ionization coefficients and equal velocities are assumed for both 
holes and electrons. These assumptions cause the E-field to be sym- 
metrical about the center of the depletion layer of a PIN diode. This 
in turn simplifies the analysis without changing the essential charac- 
teristics of TRAPATT oscillations. These assumptions also made it 
possible to develop a time-domain computer program for TRAPATT oscilla- 
tions that is two or three orders of magnitude faster than more exact 
computer programs. Voltage waveforms for the diode were computed for 
several different current waveforms using this program. Typical wave- 
forms and computer powers and impedances are given in Section 2. 
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Section 3 contains a description of a second type of high- 
efficiency mode. A trapped plasma is generated in this mode as in the 
ordinary TRAPATT mode. However, in this case, the voltage is driven 
into avalanche by the diode current, which then drops to zero (or to a 
small value) while the high field generates sufficient charge to drop 
the voltage and increase the current again. A semianalytical solution 
for the diode voltage waveform is derived for a rectangular current 
waveform. No specific calculations of efficiency or required circuit 
impedances have been made for this mode. However, it is believed that 
many experimental observations of high efficiency oscillation may be 
representative of this mode of operation. 

In Section 4 we discuss a third type of high-efficiency oscillation 
in a PIN avalanche diode. The mode is characterized by very large drive 
current densities and corresponding avalanching over a major portion of 
a half cycle of oscillation. In the modes described in Section 2 and 
3, it is assumed that the drive current does not exceed a maximum value 
so that an avalanche does not occur more than once during a cycle. In 
Section 4 it is shown that if the current is significantly larger than 
the maximum current mentioned after the initial trapped plasma is par- 
tially removed, the diode will break down again. Within a small 
fraction of a cycle the diode may reach a quasi-static state of continu- 
ous breakdown at the edges of the depletion layer while the dense trapped 
plasma in the center is being removed. When the plasma is removed, the 
current decreases and the voltage increases for the second half of the 
cycle. Preliminary calculations of this mode show that the output 
power increases proportional to the diode current with a small decrease 
in efficiency compared to the ordinary TRAPATT mode. 

Section 5 does not necessarily discuss a different mode of operation, 
but contains another potentially useful analytical analysis for 
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moderately high-efficiency oscillations in a PIN diode. By assuming 
a step function for the ionization coefficient and saturated velocities 
for the charge carriers, it is possible to divide the depletion layer 
into avalanching regions and nonavalanching regions. Correspondingly, 
by an appropriate choice of dependent variable, the telegraphers 
equation and the wave equation respectively describe the device dynamics. 
An analytical solution for one admissible pattern in the time-space 
plane of the junction is given in Section 5. It is shown that the 
efficiency can be optimized by an appropriate choice of the current 
waveform in certain portions of a cycle. 

A review of the experimental results is contained in Section 6. 

The several different circuits used are described and their relative 
merits are discussed. 
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2. TRAPATT OSCILLATIONS IN A PIN DIODE 


The anomalous mode of oscillation in avalanche diodes was first 

1 ♦ 

reported by Prager et al . This mode was characterized by high- 

efficiency (60 percent) oscillations at frequencies considerably below 

the diode avalanche frequency. Detailed computer simulations of a high 

2 

efficiency mode made by Johnston et al . showed that during part of a 
cycle of oscillation, a plasma was created in the diode depletion layer 
and was "trapped" by the low electric field. This mode was called 
TRAPATT, for Trapped Plasma Avalanche Triggered Transit. 

+ + 

Approximate analytical solutions for the TRAPATT mode in a P NN 

3 4 

diode have been developed. 9 These analyses show that a high-field 
avalanche "zone" propagates through the diode and fills the depletion 
layer with a dense plasma of electrons and holes which become trapped 
in the low-field region behind the "zone." In this section we present 
an approximate solution for the voltage waveform in a PIN diode that is 

t 

oscillating in the TRAPATT mode. It is shown that the traveling 
avalanche "zone" is not essential to generate a trapped plasma. We 
also show that the boundary conditions confine the trapped plasma to 
the central portion of the depletion region. 


References are listed at the end of the report. 

It has been shown by Muller and Guckel that the dc negative resistance 
in a PIN diode gives rise to an instability that tends to lead to the 
formation of current filaments. 5 The oscillations discussed in this 
section result from a dynamic negative resistance that can be con- 
siderably larger than the dc negative resistance and the current 
filaments should not form. 
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An approximate solution for the voltage of a PIN diode operating 
in the TRAPATT mode is derived assuming a square-wave current. A square- 
wave current was chosen for mathematical convenience although it leads 
to high-efficiency operation. Expressions for the diode voltage are 
developed for each of seven separate portions of a cycle as shown in 
Figure 2-1. The voltage waveform is very similar to that for a P + NN + 



TIME 

TA-8327-1 

FIGURE 2-1 TYPICAL VOLTAGE WAVEFORM FOR THE TRAPATT 
MODE OF AN AVALANCHE DIODE ASSUMING A 
SQUARE-WAVE CURRENT DRIVE 


diode. At point A, the E-field is uniform throughout the sample, and 
its magnitude is large but less than the value required for avalanche 
breakdown. At this instant, the diode current is turned on. Since the 
only charge carriers present are those due to thermal generation, the 
diode initially charges up like a linear capacitor, driving the magni- 
tude of the E-field above the value required for avalanche. When a 
sufficient number of carriers are generated, the particle current exceeds 
the external current, and the E-field is depressed throughout the 
depletion region, causing the voltage to decrease. This portion of the 
cycle is shown by the curve from point B to point C. During this time 
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interval^ the E-field is sufficiently large for avalanche to continue, 
and a dense plasma of electrons and holes is created. As some of the 
electrons and holes drift out of the ends of the depletion layer, the 
field is further depressed and "traps" the remaining plasma. The vol- 
tage decreases to point D. A long time is required to remove the plasma 
because the total plasma charge is large compared to the charge per unit 
time in the external current. 

At point E the plasma is removed, but a residual charge of elec- 
trons remains in one side of the depletion layer and a residual charge 
of holes in the other side. As the residual charge is removed, the 
voltage increases from point E to point F. At point F, all the charge 
that was generated internally has been removed. This charge must be 
greater than or equal to that supplied by the external current or the 
voltage will exceed that at point A. This places a limit on the 
allowable magnitude of the diode current to prevent a second avalanche 
within the same period. From point F to point G, the diode charges up 
again like a fixed capacitor. At point G, the diode current goes to 
zero for half a period, and the voltage remains constant at V until 

A 

the current comes back on and the cycle repeats. The plasma extraction, 
residual charge extraction, and final charge-up phases of the cycle are 
as described by others except that N - N = 0 in a PIN diode. 

D n 

In Sections 2-A through 2-G, approximate expressions for the vol- 
tage are derived for each portion of the waveform. These results are 
used to compute the diode power, impedance, and conversion efficiency 
at the fundamental and at several harmonics as a function of the diode 
saturation current and the magnitude of the drive current. The results 
are only indicative of the diode operation rather than exact calcula- 
tions of the diode operating characteristics. Equal ionization coeffi- 
cients for holes and electrons and the same functional dependence on E 
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for electron and hole velocities were assumed. The essential features 
of the TRAPATT mode are not obscured by these assumptions. 

A time-domain analysis computer program has also been written to 
model TRAPATT oscillations in a PIN avalanche diode when the current is 
not a square wave. This program has a running time two or three orders 
of magnitude less than more exact time-domain computer analyses but 
exhibits all the essential diode terminal characteristics. The results 
of the approximate model are compared with those of a more exact computer 
model over portions of a TRAPATT cycle, and it is shown that the approx- 
imate model is relatively accurate. The computer model is used to study 
several different diode current waveforms. Typical results for diode 
power, dc-to-RF conversion efficiency, and required circuit impedances 
are presented for different current waveforms. The obvious advantage 
of this simplified model is that parametric studies can be made 
relatively economically. 


A. Basic Equations 

The electric field, E, the electron density, n, and the hole 
density, p, in the depletion region of a P-N junction are related by 
the continuity equations and Poisson's equation. In normalized units 

Jjc 

these equations are written as follows: 

d n d J n 

" aT + “ v [n + pI <2 - 1) - 


. a j 

O P P r , 

— = - + Q(v [n + p] 

o t a x 


(2-2) 


* 

The normalization is given in Appendix A. 
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— = N -N + p - n 
9 x D A 


(2-3) 


where J = vn, J = vp, and a(E) is the ionization coefficient and v is 
n P 

the particle velocity. The total external current density, J , can be 
derived from Eqs. (2-1), (2-2), and (2-3): 

— = - (J + J )' + J . (2-4) 

9 1 n p 

The boundary conditions at the edges of the depletion layer are 

J (0) = J (x)| = J , (2-5a) 

p p 1 x=0 s 

J (1) = J (x) I = J . (2-5b) 

n n n=l s 

In a PIN diode, N - N is zero, and in the absence of internal charge, 

D A 

the E-field is uniform. When the charge density is small or moderate 
and the average value of the field is high, the E-field is still 
essentially uniform. When the depletion layer is filled with a plasma, 
the E-field is uniform throughout the plasma except for a finite 
gradient near the edges. Therefore, over certain portions of a TRAPATT 
cycle it is useful to write the equations for a PIN diode under the 
assumption that the E-field is uniform. To account for the removal of 
charge by the carriers drifting out of the depletion region, we will 
introduce a charge relaxation. 

Equations (2-1) and (2-2) can be combined to give 

— = — (J - J ) + 2 avQ , (2-6) 

9 t 9 x p n 
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where Q = n + p and is equal to the total charge density. Integrating 
with respect to x, and taking a and v outside the integral since E is 
effectively uniform, we obtain 


dQ , v , 1 

— = - (J - J ) + 2 avQ , 

at ■ p n ' 


(2-7) 


where Q 


Qdx. If the E-field were exactly uniform, J - J would 
0 p n 


be zero, but due to the boundary conditions, J (0) = J (l) = J . We 

p n s 

approximate J (1) and J (0) by avQ, where a is some proportionality 
* P n 

constant : 


dQ 

— = 2v(of - a)Q + 2J 
dt s 


(2-8) 


The charge removed by drift from the depletion layer has been modeled by 

a relaxation of the total charge. Obviously, this is not strictly 

correct because all the charges could drift out and J and J would be 

n p 


zero 


or J . 
s 


If we neglect J , 


the solution of Eq. 


(2-8) is 


Q q exp 


) 2 p t 

Po 


v(a - a) dt! 


(2-9) 


where Q q is the charge density at t = 0. Integrating Eq. (2-4) with 
respect to x and assuming a uniform E-field yields: 


It can be shown on a small-signal basis that a more exact approximation 
for Jp(l) and J n (0) would be avQ(t - t), where t corresponds to a drift 
time. For the TRAPATT mode, however, the former approximation is 
adequate and somewhat simpler. 
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— = - vQ + J 
dt 


(2-10) 


Note that for a uniform field, the normalized terminal voltage is equal 
to the magnitude of the E-field. Most of the characteristics of the 
TRAPATT mode for a PIN diode are exhibited by the simple model of an 
avalanche diode given by Eqs. (2-8) and (2-10). This model is inade- 
quate only during that portion of a cycle when the trapped plasma is 
being extracted. Since the model assumes that the plasma is removed 
uniformly across the sample, it does not account for the electric field 
being larger in those regions of the depletion layer where the plasma 
has been removed. In Appendix B we show how to correct for the plasma 
removed from the ends of the depletion layer, and a computer program 
is described that numerically integrates Eqs. (2-8) and (2-10). Results 
from this computer program have shown that the proportionality constant, 
a, is approximately unity in order for all the charge generated during 
avalanche to equal the integral of the external current over one cycle. 

B. Initial Charge-Up Phase 

In Figure 2-1, the E-field that corresponds to the voltage at point 
A is uniform and just below avalanche. Initially the charge density is 
small, so Eqs. (2-8) and (2-10) can be used to describe the field and 
the voltage from point A to point B. For an applied constant current 
density J, the E-field initially increases linearly with time: 

E = E + Jt . (2-11) 

A 

It continues to increase approximately linearly until vQ equals J, and 
the voltage has reached its maximum value. From Eq. (2-9) we can 
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determine the time, t , when the voltage reaches its maximum value, if 

m 

we use an exponential function for a : 

XE 

ot = a e . (2-12) 

o 


The constants a and X are chosen by fitting the exponential function 



Figure 2-2 shows a fit to the data for electrons in silicon. 



1.0 1.2 1 4 1.6 1.8 2.0 2.2 2.4 


NORMALIZED E-FIELD 

TA— 8327— 4 

FIGURE 2-2 COMPARISON OF EXPONENTIAL APPROXIMATION 

OF THE IONIZATION COEFFICIENT FOR ELECTRONS 
IN SILICON WITH A MORE EXACT FORM 


Substituting Eq. (2-11) into Eq. (2-12) and using the result in 
Eq. (2-9 ) } we obtain 


Q = Q 


s 



(2-13) 
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where Q = J /v and a = a(E ). According to our normalization, v = 
s s s A A s 

1; we retain v in our equations for completeness and clarity. We have 
s 

set the constant, a, equal to unity for the reason discussed earlier. 

We obtain a transcendental equation for t by setting Q = J/v : 

m s 


t = £n |l + — \t + 1/ (2v X,n(J/J . (2-14) 

mXJ / m s S J) 


An iterative procedure is used to evaluate t . Initially, t is set 

m m 

equal to zero in the argument of the logarithm, and a first approxi- 


mation for t is obtained. This value is used to obtain a better 
m 

approximation. Two or three iterations are sufficient if the ratio 


J/J is large, 
s 


The voltage from point A to point B is approximately 


V(t) — V + J t , 0 £ t £ t . (2-15) 

A m 


C. Plasma Generation 

During the time from point B to point C in Figure 2-1, a large 
concentration of holes and electrons is generated in the depletion 
layer. An accurate calculation of the total charge generated and the 
corresponding spatial distribution of the E-field requires a computer 
program that integrates the partial differential equations (2-1), (2-2), 
and (2-3). We have written such a program, and it will be referred to 
as UNSAT. Examination of some computer simulations shows that during 
the short time from point B to point C, the E-field is effectively 
uniform across the sample except near the boundaries. Therefore, 

Eqs. (2-9) and (2-10) can also be used during this portion of the 
TRAPATT cycle. 
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To obtain an analytical solution for Eqs. (2-9) and (2-10) over the 

time interval, t (from point B to point C) , we assume that a is equal 
g 

to a constant = a(E ) for values of E greater than some critical 
m max 

field value, E , and is zero for E < E . This approximation is similar 
c c ++34 

to that used in the analytical solutions of a P NN diode. y 

The solution to Eq. (2-8) for constant a and v = v is 

s 


2v to -l)t 

Q(t) = Q e S m 
m 


(2-16) 


where Q = J/v is the charge density at maximum E-field. Substituting 
m s 

Eq. (2-16) into Eq. (2-10), we obtain for E: 


E = E + J(t-t ) - — 

m m 2v to -1) 

s m 


2v to -l)(t-t ) 
s m m 

e - 1 


t £ t <£ t + t . (2-17) 
m g m 


When E = E , a is zero, and the avalanche process stops. The time, t , 
c g 

required to go from point B to point C can be found from Eq. (2-17) by 

setting E = E , and t = t + t : 
c m g 


t 

g 


2v to -l) 
s m 


2v to -1) f 

In { 1 + — - — E - E + Jt 

J L m c g^ 


(2-18) 


One can solve for t numerically by iteration, as explained earlier 

g 


for t . 
m 


It we use the right-hand side of Eq. (2-18) in Eq. (2-16) we 

obtain an expression for the total charge, Q , in the plasma filling 

P 

the depletion layer: 
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(2-19) 


Q = Q + 2 (a -1) E - E + Jt 
p m m L m c , g_ 


The total number of holes or electrons generated in the time t 

is G: 

t + t 

G =Jl g m vaQdt . (2-20) 

t 

m 


Equation (2-20) can be written as 


G = oi |E - E + Jt . (2-21) 

mL m c gJ 


An additional negative or positive charge of 1/2 Q was generated while 

m 

the E-field was increasing to its maximum value. 

The importance of G is that it must equal the integral of the' 
external current over one period. This condition is used as a check 
that the times calculated in the approximate solution all add up to the 
correct period for a square wave. This condition is even more useful 
for the non-square-wave current waveforms discussed later. 


D. Field Depression and Partial Plasma Extraction 

We can rewrite Eq. (2-19) in the following form: 

Q = 2G - 2|E - E + Jt | . (2-22) 

p L m c gJ 

This equation shows that the plasma charge density throughout the 
depletion layer is equal to the total charge generated less a quantity 
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equal to the charge that has drifted out of the depletion layer. It is 
the drifting charges that actually depress the E-field and stop the 
avalanche process. The amount of charge that drifts out is equal to 
that amount required at the edges of the depletion region to drop the 
field from E to E plus the charge in the external current during; 
generation. 

Beginning at point C in the cycle, we can no longer use the assump- 
tion of a uniform E-field. The boundary conditions do not permit a 
uniform E-field throughout the depletion layer in the presence of large 
charge densities. At the left boundary the hole current is equal to 
the saturation current, while the electron current has a large finite 
value and vice versa at the right-hand boundary. On the other hand, the 
field inside the plasma region is uniform and equal to some small value 
determined by the plasma velocity. These two conditions are satisfied 
simultaneously by the charge carriers in the plasma near the boundaries 
leaving the depletion layer. A residual charge of one carrier type 
remains in each end of the depletion layer sufficient to satisfy Poisson 1 s 
equation and the total current equation. 

Initially, when the E-field is still high the carriers in the plasma 

move at saturated velocity. As the field decreases, the carrier velocity 

decreases until it reaches the plasma velocity, v , given by 

P 

V P = J/Q P * (2-23) 

where Q is given by Eq. (2-22). In the plasma the total change in the 

ir 

E-field is E - E , where E is the plasma field. If t is the time 
c P P d 

required for the E-field to change by this amount, then we have, 

E - E + Jt = Total charge removed during t . (2-24) 
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We assume that all the charge is removed from the ends of the depletion 

layer within a depth, x , from each boundary. The number of holes per 

d 

unit distance removed near the left boundary is 1/2 Q ; the number of 

P 

electrons per unit distance removed near the same boundary is 

(1/2 Q - n ), where n is the residual electron density; thus from 
p r r 

Eq. (2-24) we have: 


x 

E - E +Jt = f d (Q - n )dx . (2-25) 

c p d o n r 


The residual charge would be equal to 



( 2 - 26 ) 


which is a function of both time and position. Only the final value of 

the residual charge density and the value of x when the field in the 

d 

plasma has dropped to E are of primary interest because this phase of 

P 

the TRAPATT cycle happens quickly. These quantities determine the 
minimum voltage at point D and the amount of charge that must be removed 
at the plasma velocity. We assume that at point D, n^ is uniform in 
the region from the boundary to the plasma. Then the E-field has the 
following distribution at the left edge of the depletion layer: 


E=E + n(x -x), O^x^x 
p r d d 


(2-27) 


A similar expression can be written for the E-field at the right boundary. 

Since x = v t, the partial derivative with respect to time of 
d p 

Eq. (2-27) when substituted into Eq. (2-26) gives 
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( 2 - 28 ) 


n = J/ (v + v ) 
r p 


The velocity, v, depends upon the E-field, which in turn is a function 
of n^_. We will assume v = as a first approximation; then Eq. (2-26) 
becomes 


E - E + Jt 
c p d 


Q 

P 


1 


(v + v ) 
s p J 


X 


d 


(2-29) 


An estimate of t , can be obtained from the fact that the holes 
d 

traverse the distance x with some average velocity, v . The initial 

d a 

velocity is v and the final velocity at point D is v ; we use a simnle 
s P 

average as a first approximation: 


2x 

d 


d v + v 
s p 


and 


E - E 

X = S E 

d , 2i_ 

p (v + V ) 
s p 


(2-30) 


(2-31) 


We have not developed an expression for the voltage during the transition 
from point C to point D. This time is short and the voltage variation 
large; we assume simply a linear relationship: 


V(t) = V 

c 


t-(t + t ) 

„ V mg 

(v “ V ) , t + t ^ t £ t + t + t 

C P t mg mgd 


(2-32) 


where 


2 

V = E + n x 
p p r d 
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This approximation should not affect any low-frequency components of the 
voltage waveform. 


E ■ Plasma Extraction 

The plasma is extracted by the external current as the voltage moves 

from point D to point E. During this time, the E-field inside the plasma 

is time-invariant and uniform. The edges of the plasma move toward the 

center of the depletion layer with the velocity v and leave a residual 

P 

charge of electrons on one side of the. depletion layer and of holes on 

the other side. The plasma exists until the right- and left-hand edges 

of the plasma meet at the center of the depletion layer. The time, t , 

P 

required to remove the plasma is 


Q (0.5 - x ) 

t S- . (2-33) 

P j 


As the plasma is extracted, the electric field between the left edge of 

the depletion layer and the left edge of the plasma, x , is given by 

P 


E = E • + n (x -x) , 0 £ x £ x , (2-34) 

P r p p 


and the voltage across the device is 


2 

V = E + n x 
P r p 


(2-35) 


The distance, x = v t, increases linearly with time until it equals 0.5. 
P P 

If t = t + t + t and t = t + (0.5 - x )/v , then the voltage is 
1 m g d 2 1 dp 
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2 

Jx 

d 


(2-36) 


V(t) 


E + 
P 


v 

s 


+ V 

p 


+ 


2 

Jv 

p 


V + v 
s p 




t < t £ t 
1 2 


> 


while the plasma is being removed. 

It was assumed for the voltage waveform shown in Figure 2-1 that 

avalanche occurred only once during a cycle. To prevent a second 

avalanche, the current J must not exceed a maximum value, J . The 

max 

maximum current can be obtained from a consideration of Eq . (2-34). At 

x = 0, E cannot exceed E , the value of the field at the beginning of 

n 

the cycle, or 


E 


A 


E > 
P 


J x 

P 


(v + V ) 
s p 


(2-37) 


The field reaches its peak value when x =0.5; thus, 

P 


J £ 2 (E - E ) (v + v ) , (2-38) 

A p s p 


4 

which agrees with DeLoach and Scharfetter. 


F . Residual Charge Extraction 

When the plasma is completely removed at point E, a residual charge 
of electrons exists on the left-hand side of the depletion layer and a 
similar residual charge of holes exists on the right half of the deple- 
tion layer. The time required to remove the residual charge is of the 


This second avalanche is the same as the "premature avalanche" dis- 
cussed by Clorfeine et al . 3 
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order of half a transit time. As the residual charge is removed, the 

field in the center increases. From Eq. (2-34) and the fact that x = 

P 

0.5 when the voltage is at point E, we can write for the field on the 
left side of the depletion region: 


E(x) = E + n (0.5 - x) , 0 ^ x ^ x , (2-39) 

p r r 


where x is the extent of the residual charge from the left boundary, 
r 

The value of x decreases linearly with time: 
r 


x = 0.5 - v (t-t ) , t £ t £ t , (2-40) 

r s 2 9 2 3 9 


where 


t 

3 


t + l/2v . The corresponding terminal voltage is 
2 s 


J(1 + 2v (t-t )) 

V(t) = E + “7 r , t £ t £ t . (2-41) 

p 4 (v + v ) 2 3 

s p 


The voltage increases linearly with a slope approximately equal to one- 
half that for a charge-free region. 

G. Final Charge-Up Phase 

At point F all the charge is removed and the field is uniform across 
the sample again. Since the field is below avalanche breakdown, the 
diode charges up like a linear capacitor. The voltage is 

V(t) = E + — — + J(t-t ) , t s t £ T/2 (2-42) 

p 2 (v + v ) 3 3 

s p 
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The time required for this final charge-up phase depends on the value 
of the voltage at point F. When J equals the maximum allowable current, 
the voltage reaches just as the last of the residual charge leaves 
the depletion layer and no additional time is required. For smaller 
currents, additional external charge is necessary to charge the diode 

back up to V . 

F A 

H * Calculations 

The semianalytical solution has been used to calculate the perfor- 
mance of a PIN diode with a 10-micron-wide depletion layer. The 
normalized low-field velocity, v, was assumed to be ten times the 
electric field. The expression of Scharfetter and Gummel for the 
electron ionization coefficient in silicon was assumed to correspond to 
the experimental data for the ionization coefficient. 6 

Figure 2-3 shows the normalized fundamental impedance of a PIN 
diode operating in a TRAPATT mode for three levels of a square-wave drive 
current. The different levels of the saturation current correspond to 
different thermal generation rates of carriers in the diode; the lower 
the level of the saturation current, the smaller the thermal generation. 
It is evident from Figure 2-3 that the higher the level of the drive 
current, the smaller the variation in the fundamental impedance over 
several orders of magnitude of the saturation current. 

Tables 2-1, 2-2, and 2-3 show the normalized period, T, power, P , 
and impedance, Z , corresponding to the results shown in Figure 2-3. 

The fundamental frequency for each data point is different; in fact, 
for a given drive current, the frequency varies over an octave, with the 
range of saturation currents shown. Also, the diode impedances at 
higher harmonics change appreciably. These results for a square-wave 
current indicate that at the higher harmonics, the diode furnishes 
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FIGURE 2-3 NORMALIZED FUNDAMENTAL IMPEDANCE OF PIN 
DIODE OPERATING IN A TRAPATT MODE WITH A 
SQUARE-WAVE CURRENT DRIVE 


power for small saturation currents and high drive levels, but that 
ac power must be furnished to it for large saturation currents and low 
drive currents. The choice of a square-wave drive current was arbitrary 
and made to simplify the analysis in the analytical solution. In 
practice, the diode would always be placed in a passive circuit, and 
the diode impedances would change accordingly. 

The results shown in Figure 2-3 and Tables 2-1, 2-2, and 2-3 are 
based upon a particular choice of parameter values used in the semi- 
analytical solution. 

The parameter, R, is the ratio of the maximum value of the ioni- 
zation coefficient, a , to that value of the ionization coefficient, a , 

m c 

below which we assume that avalanche stops and the ionization coefficient 
drops to zero. By comparing the results of a more exact computer 
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Table 2-1 


NORMALIZED PERIOD, T, POWERS, P., AND 

IMPEDANCES, Z , FOR A PIN AVALANCHE DIODE OPERATING IN 

THE TRAPATT MODE, WITH A SQUARE-WAVE CURRENT DRIVE, 

AS COMPUTED BY THE ANALYTICAL SOLUTION 

(J = 1.5, E = 1.0, R = 2.5) 

A 


Normalized 

Parameter 

-9 

J = 2 X 10 
s 

-7 

J = 2 X 10 
s 

-5 

J = 2 X 10 
s 

-3 

J = 2 X 10 
s 

Period 

25.5 

20.0 

15.4 

10.0 

P 

0 

0.473 

0.483 

0.500 

0.535 

P 

1 

-0.266 

-0.263 

-0.255 

-0.234 

P 

3 

-0.019 

-0.015 

-0.009 

+ 0. 005 

P 

5 

-0.002 

-0.001 

+0.006 

+0.011 

P 

7 

+0.002 

-0.004 

+0.006 

+0.005 

Z 

0 

0.840 

0.858 

0.888 

0. 951 

Z 

1 

-0.584 + j 0.058 

-0.576 + j 0.080 

-0. 560 + j 0.112 

-0.513 + j 0. 168 

Z 

3 

-0.374 + j 0.306 

-0.304 + j 0.336 

-0.176 + j 0.360 

+0.105 + j 0.316 

Z 

5 

-0.083 + j 0.406 

+0.070 + j 0.383 

+0.310 + j 0.275 

+0.630 - j 0.120 

Z 

7 

+0.252 + j 0.346 

+0.447 + j 0.206 

+0.657 - J 0.110 

+0. 534 - j 0.707 



Table 2-2 


NORMALIZED PERIOD, T, POWERS, P., AND 

IMPEDANCES, Z , FOR A PIN AVALANCHE DIODE OPERATING IN THE TRAP ATT 
i 

MODE, WITH A SQUARE -WAVE CURRENT DRIVE, AS COMPUTED 

BY THE ANALYTICAL SOLUTION 

(J = 1.0, E = 1.0, R = 2.4) 

A 


Normalized 

Parameter 

-9 

J = 2 X 10 
s 

-7 

J = 2 X 10 

s 

-5 

J = 2 X 10 
s 

-3 

J = 2 X 10 
s 

Period 

23 . 9 

19.6 

15.0 

10.6 . 

P 

0.323 

0.333 

0.349 

0.384 

0 





P 

-0.177 

-0.172 

-0.162 

-0.134 

1 





P 3 

-0.008 

-0.004 

+0 . 003 

+0.014 

P 

+0.003 

+0.005 

+0.007 

+0.006 

5 





P 

+0.004 

+0.004 

+0.003 

-0 . 002 

7 





z o 

1.291 

1.331 

1.398 

1.53 

z 

1 

-0.875 + j 0.182 

-0.849 + j 0.218 

-0.797 + j 0.271 

-0.659 + j 0.353 

z 

3 

-0.365 + j 0.537 

-0.185 + j 0.54 

+0.118 + j 0.471 

+0 . 636 + j 0 . 109 

Z 

5 

+0.327 + j 0.454 

+0.611 + j 0.242 

+0.910 - j 0.219 

+0.710 - j 1.052 

z 7 

+0.879 - j 0.060 

+1.000 - j 0.525 

+0.755 - j 0.193 

-0.481 - j 0.973 
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Table 2-3 


NORMALIZED PERIOD, T, POWERS, P., AND 
IMPEDANCES, Z , FOR A PIN AVALANCHE DIODE OPERATING IN THE TRAPATT 
MODE, WITH A SQUARE-WAVE CURRENT DRIVE, AS COMPUTED 
BY THE ANALYTICAL SOLUTION 
(J = 0.8, E = 1.0, R = 2.4) 

n 


Normalized 

Parameter 

-9 

J = 2 X 10 
s 

-7 

J = 2 X 10 
s 

-5 

J = 2 X 10 
s 

-3 

J = 2 X 10 
s 

Period 

25.0 

20.5 

16.0 

11.2 

P 

0 

0.259 

0.266 

0.278 

0.303 

P 

1 

-0. 142 

-0. 138 

-0.130 

-0.111 

P 

3 

-0.006 

-0.003 

+0.002 

+0.011 

P 

5 

+0.003 

+0.004 

+0.006 

+0.005 

P 

7 

+0.003 

+0.001 

+0.002 

-0.002 

Z 

0 

1.617 

1.664 

1.740 

1.893 

z 

1 

-1.096 + j 0.238 

-1.064 + j 0.276 

-1.005 + j 0.330 

-0.852 + j 0.412 

z 

3 

-0.404 + j 0.657 

-0.188 + j 0.643 

+0.167 + j 0.546 

+0.765 + j 0.138 

Z 

5 

+0.505 + j 0.467 

+0.821 + j 0.190 

+1.136 - j 0.348 

+0.914 - j 1.271 

Z 

7 

+1.147 - j 0.282 

+1.220 - j 0.822 

+0.861 - j 1.467 

-0.587 - j 1.226 




simulation (UNSAT) with those of the analytical solution, over the 
initial charge-up and generation phases of a TRAPATT cycle, we found 
that the value of R differs for different drive currents. Figure 2-4(a), 
(b), and (c)" compare the amount of total charge generation as a function 
of the saturation current as calculated by the computer simulation 
(UNSAT) and by the analytical solution. Note that for J = 0.9, R is 
2.4; whereas for J = 1.0 and 1.5, R is 2.4 and 2.5, respectively. Since 
the diode fundamental period of oscillation for a square-wave drive 
current is equal to twice the generation divided by the current magni- 
tude, Figure 2-4 can be used to determine the frequency. 

We have also made a comparison of the diode voltage as a function 
of time as computed by the approximate analytical solution and by the 
exact computer simulation (UNSAT) over the initial charge-up and charge 
generation phases of the TRAPATT cycle. These results are shown in 
Figure 2-5(a) and (b) for two different drive levels. The results from 
the approximate analytical solution are quite good. A comparison was 
not made over a longer portion of the cycle because of the excessive 
running times required for the exact computer simulation. 

I . Nonsquare-Wave Current Waveforms 

Several other current waveforms were investigated using the computer 
program modified SMPAL described in Appendix B. For example, current 
waveforms composed of the first few nonzero harmonics of a square-wave 
were used. The relative magnitudes and phases of the harmonics were 
the same as for a square wave. Figures 2-6 and 2-7 show the normalized 
voltage waveforms for currents composed of the first and third harmonics 
and of the first, third, and fifth harmonics, respectively. 

Table 2-4 lists the powers and diode impedances corresponding to 
Figures 2-6 and 2-7. Comparison of Table 2-4 with Tables 2-1 and 2-2 


27 




NORMALIZED SATURATION CURRENT — J $ 
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FIGURE 2-4 TOTAL NORMALIZED CHARGE GENERATION DURING ONE 
CYCLE OF THE TRAPATT MODE IN A PIN DIODE 
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ANALYTICAL SOLUTION 





FIGURE 2-6 NORMALIZED VOLTAGE AND CURRENT WAVEFORMS FOR 
THE TRAPATT MODE WITH A CURRENT USING THE FIRST 
AND THIRD HARMONICS OF A SQUARE WAVE 


shows that for the same and J, the fundamental period is about the 
same for the nonsquare-wave waveform. This is a consequence of the 
fact that the currents in Figures 2-6 and 2-7 do charge the diode up 
about as rapidly as does a squarewave current. For the assumed 
operating conditions in Figures 2-6 and 2-7, power must be fed to the 
diode at the third and fifth harmonics. 
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Table 2-4 


NORMALIZED PERIOD, T, POWERS, P^, AND 

IMPEDANCES, Z , FOR A PIN AVALANCHE DIODE FOR A NORMALIZED 

DRIVE CURRENT USING UP TO THE SEVENTH HARMONICS OF A 

SQUARE WAVE, AS COMPUTED BY MODIFIED SMPAL 

(V = 1.0, J = 2 X lO -9 ) 

A s 

J = A [0. 5 + 0.6366 sin 9 + 0.2122 sin 30 + 0.1273 sin 50 + 0.0909 sin 70] 


Normalized 

Parameter 

Three 
Harmonics 
A = 1.0 

Five 

Harmonics 
A = 1.5 

Seven 
Harmonics 
A = 1.5 

Period 

24.57 

24.90 

26.4 

p o 

0.26 

0.383 

0.41 

p 

-0.10 

-0.15 

-0.16 

1 




p 

3 

+0.004 

-0.0002 

-0.006 

P 


+0. 007 

+0.004 

5 




P 

7 

— 

— 

+0.005 

Z 

0 

1.04 

0.68 

1 

0.73 

z 

1 

-0.51 + j 0.24 

-0.32 + j 0.08 

-0.35 + j 0.04 

z 

3 

+0.17 + j 0.66 

-0.003 + j 0.40 

-0.11 + j 0.44 

Z 

5 

— 

+0.40 + j 0.43 

+0.20 + j 0.44 

Z 

7 

— 

— 

+0.50 + j 0.32 
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FIGURE 2-7 NORMALIZED VOLTAGE AND CURRENT WAVEFORMS FOR 
THE TRAPATT MODE WITH A CURRENT USING THE FIRST, 
THIRD, AND FIFTH HARMONICS OF A SQUARE WAVE 


The efficiency has also decreased appreciably from that for a 
square wave. This is primarily a consequence of the voltage not being 
as high in the second half of the cycle as in the case of a square 
wave. The voltage is less in this region because the charge necessary 
to bring the voltage back up to the original starting value is supplied 
by the current at the end of the cycle. 
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The efficiency can be improved if the entire current waveform is 
shifted to the right and all other parameters remain the same. This is 
shown in Figure 2-8 for the case of a current composed of the. first and 
third harmonics of a square wave. Note that the normalized voltage 
during the last half of the cycle now oscillates around unity instead of 
about 0.8. The dc-to-RF conversion efficiency in the fundamental is 50 
percent, and the real part of the diode impedance at the third harmonic 
is negative. Similar shifts have been made for current waveforms using 
up to the fifth and up to the seventh harmonics of a square wave. These 



6 — degrees 


TA-65 1524-4 


FIGURE 2-8 NORMALIZED VOLTAGE AND CURRENT WAVEFORMS FOR 
THE TRAPATT MODE WITH A CURRENT USING THE FIRST 
AND THIRD HARMONICS OF A SQUARE WAVE SHIFTED 
15 DEGREES 
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results are tabulated in Tables 2-5 and 2-6. The efficiency has increased 
over the corresponding results for the unshifted curves. The impedances 
at the higher harmonics either are effectively reactive or have a small 
negative real part except for the case of the seventh harmonic. 

The shifting of the current waveform is equivalent to changing the 

impedance of a passive circuit. These results indicate, therefore, that 

by proper tuning of the circuit, 50-percent efficiency should be possible 

with just the first few harmonics of a square-wave current. This is in 

+ + 3 

agreement with the results for a P NN diode. 

If the initial voltage, V^, at time t = 0, were increased, the 

effect would be similar to shifting the current waveform to the right. 

The two effects would be exactly equivalent if the initial charge were 

also increased above the thermal generation value to account for the 

amount of generation that occurs from V = 1,0 to the increased value of 

A 

voltage. 

Other current waveforms were tried where for the assumed operating- 
conditions a trapped plasma occurred, but the diode exhibited no negative 
resistance. In each of these cases, too much of the charge required 
from the external current to neutralize the avalanche charge occurred 
during the second half of a period. Thus, the voltage remained low 
after avalanche until the very last portion of the cycle. The higher 
harmonics are very important in shaping the current waveform so that it 
is nearly zero over the entire last half of a cycle. 

The half-wave rectified sine wave was considered as a drive current 
waveform because the current is zero over the second half of a cycle. 

Also, there are no discontinuities in this waveform, so the higher 
harmonics are less important than in a square wave. In fact the first 
and second harmonics are the most significant. A normalized current 
density of the following form was used: 
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Table 2-5 


NORMALIZED PERIOD, T, POWERS, P i , AND IMPEDANCES, Z^, FOR A PIN AVALANCHE 

DIODE OPERATING IN THE TRAPATT MODE 

j’ = A [0. 5 + 0.6366 sin (9 - 15°) + 0.2122 sin (39 - 45°)], V, = 1.0 

A 



J = 2 
s 

-9 

x 10 


Parameter 

A = 1.0 

A = 1.5 

A = 2.0 

Period 

22.5 

20.28 

18.84 

P 

0 

0.35 

0.52 

0.71 

P 

1 

-0.17 

-0.26 

-0.32 

P 3 

-0.005 

-0.01 

-0.02 

z o 

1 .42 

0.93 

0.71 

z 

1 

-0.86 + j 0.19 

-0.57 + j 0.06 

-0.40 + j 0.01 

z 

3 

-0.22 + j 0.63 

-0.30 + j 0.39 

-0.24 + j 0.28 





Table 2-6 


NORMALIZED PERIOD, T, POWERS, P^ AND IMPEDANCES, Z., FOR A PIN 

AVALANCHE DIODE FOR A NORMALIZED DRIVE CURRENT USING UP TO FIFTH AND 

UP TO SEVENTH HARMONICS OF A SQUARE WAVE SHIFTED 1/36 PERIOD TO THE RIGHT 

V = 1.0, A = 1.5, J = 2 X lO -9 
A s 


Parameter 

Five Harmonics 

Seven Harmonics 

Period 

23.0 

24.70 

p o 

0.52 

0.52 

p 

1 

-0.26 

-0.26 

p 

3 

-0.015 

-0.018 

P 

5 

0.001 

-0.001 

P 

7 

— 

0.002 

z o 

0.92 

0.91 

z 

1 

-0.57 + j 0.05 

-0.58 + j 0.04 

z 

3 

-0.30 + j 0.36 

-0.36 + j 0.32 

z 

5 

0.06 + j 0.46 

-0.07 + j 0.45 

Z 

7 

— 

0.27 + j 0.43 
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(2-43) 


J(t) 



/ 2nt \ 

/ 4xct o\“ 

A 

0.318 + 0.5 sin l - ^ - 1 

+ 0.212 sin + 270 1 


Figure 2-9 shows a typical voltage waveform for a current given by 
Eq. (2-43) as computed by modified SMPAL. Three different values of the 
current drive level, A, and four different levels of the saturation 
current were considered. These results are summarized in Table 2-7. 

The diode exhibits a negative resistance at both frequencies over the 



FIGURE 2-9 NORMALIZED VOLTAGE AND CURRENT WAVEFORMS FOR 
THE TRAPATT MODE WITH A CURRENT USING THE FIRST 
AND SECOND HARMONICS OF A HALF-WAVE RECTIFIED 
SINE WAVE 
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Table 2-7 


NORMALIZED PERIOD, T, POWERS, P , AND IMPEDANCES, Z., FOR 
' i i 

A PIN AVALANCHE DIODE OPERATING IN THE TRAP ATT MODE FOR A DRIVE CURRENT USING 
THE FIRST AND SECOND HARMONICS OF A HALF-WAVE RECTIFIED SINE WAVE 


Parameter 

A = 2.0 

A = 1.5 

o 

(H 

II 

< 


J 

s 

-9 

= 2 x 10 


T 

20.64 

22.49 

25.41 

P 

o 

0.500 

0.375 

0.251 

P 

1 

-0.20 

-0.15 

-0.096 

P 2 

-0.06 

-0.05 

-0.031 

z 

o 

1.24 

1.65 

2.48 

Z 

1 

-0.406 + j 0.112 

-0.540 + j 0.174 

-0.765 + j 0.300 

z 

2 

-0.660 + j 0.226 

-0.905 + j 0.390 

-1.37 + j 0.756 


J 

s 

-7 

= 2 x 10 


T 

17.65 

19.22 

21.72 

P 

o 

0.500 

0.380 

0.258 

P 

1 

-0.187 

-0.137 

-0.084 

P 2 

-0.028 

-0. 045 

-0.029 

Z 

! o 

1.24 

1.67 

2.55 

Z 

1 

-0.375 + j 0.114 

-0.491 + j 0.186 

-0.703 + j 0.304 

Z 

2 

-0.658 + j 0.248 

-0.905 + j 0.442 

-1.33 + j 0.813 


J : 
s 

= 2 x 10 5 


T 

14.41 

15.68 

17.66 

P 

o 

0.519 

0.395 

0.268 

P 

1 

-0,169 

-0.123 

-0.075 

P 2 

-0.059 

-0. 044 

-0.028 

Z 

o 

1.28 

1.75 

2.64 

z i 

| -0.338 + j 0.118 

j -0.441 + j 0.183 

-0.602 + j 0.307 

Z 2 

-0.655 + j 0.288 

-0.869 + j 0.48 

| -1.23 + j 0.887 


J = 

s 

= 2 x 10” 3 


T 

10.74 

11.63 

13.03 

P 

o 

0.550 

0.418 

0.287 

P 

1 

-0.130 

-0.088 

-0.046 

P 2 

-0.054 

-0.036 

-0.020 

z 

o 

1.36 

1.84 

2.85 

Z 

1 

-0.259 + j 0.120 

-0.316 + j 0.187 

-0.371 + j 0.306 

z 

2 

-0.598 + j 0.36 

-0.731 + j 0.578 

-0.871 + j 1.01 
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the complete range of parameters used. The reactive component of the 

impedance at each frequency is inductive. As the saturation current 

-9 -3 

increases from 2 X 10 to 2 X 10 , the fundamental frequency increases 
by approximately a factor of two, and the efficiency decreases by about 
the same factor. 

If a sinusoidal current waveform is shifted in phase and all other 
parameters kept the same, high-efficiency oscillations are possible. We 
have considered normalized diode current waveforms of the form 

J(0) = 1 . 5 [0 . 5 + 0.55 sin (9 - cp ) ] (2-44) 

and computed the corresponding voltage, power, efficiencies, and impe- 
dances of a PIN avalanche diode using the computer program, modified 

SMPAL . Figure 2-10 shows the current and’ voltage waveforms when cp = 
o 

45 . The theoretical efficiency of this case is 44 percent. Results 

for current waveforms with other phase angles are summarized in Table 

-9 

2-8, where the saturation current was assumed to be J = 2 X 10 

s 

It is unlikely that this mode has been observed experimentally 
because it is difficult to construct a microwave circuit that allows 
current to flow at only a single frequency. 

Current waveforms of avalanche diodes operating in a high-efficiency 

7 

mode have been measured by Snapp et al . We have Fourier analyzed the 
measured current waveform for one diode that Snapp reported to be 
oscillating at 49-percent efficiency. This current was then used as 
the driving current in the computer program, modified SMPAL , to compute 
the corresponding voltage, powers, efficiencies, and impedances. 

Figures 2-11 and 2-12 show the normalized voltages and current waveforms 
where the current was composed of the first three harmonics and first 
ten harmonics, respectively. Table 2-9 summarizes the operating 
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FIGURE 2-10 NORMALIZED VOLTAGE AND CURRENT WAVEFORM IN 
A PIN AVALANCHE DIODE WHERE THE CURRENT HAS 
A dc AND A FUNDAMENTAL COMPONENT 


Table 2-8 

NORMALIZED PERIOD, T, POWERS, P^ AND IMPEDANCES, Z. 
FOR A PIN AVALANCHE DIODE OPERATING IN THE TRAPATT MODE 
FOR A SINUSOIDAL DRIVE CURRENT 



T 

P 

0 

P 

1 

z 

0 

z 

i 

0 

o 

18.98 

0.19 

+0.037 

0.35 

+0 . 11+ jO . 34 

20 

o 

16.56 

0.37 

-0.12 

0.65 

-0. 35+ jO. 33 

45 

o 

14.85 

0.49 

-0.22 

0.87 

-0. 64+ jO. 28 

65 

14.49 

0. 51 

-0.23 

0. 91 

-0. 68+ jO. 26 
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FIGURE 2-11 NORMALIZED VOLTAGE AND CURRENT WAVEFORMS FOR 
THE TRAPATT MODE IN A PIN DIODE WITH A DRIVE 
CURRENT USING THREE HARMONICS OF A WAVEFORM 
REPORTED BY SNAPP, et al. 7 


characteristics for these two cases. The normalized thermal generation 
-9 

was 2 X 10 . It is evident from these results that only the first 

three harmonics of the current are significant. From a circuit point 
of view this is encouraging because it is difficult to control more 
than a few harmonics of current or voltage. 

The theoretical results for the current waveforms in Figures 2-11 

and 2-12 predict efficiencies comparable to that reported by Snapp 
7 

et al. The period of oscillation is about double. There are several 
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2.5 



FIGURE 2-12 NORMALIZED VOLTAGE AND CURRENT WAVEFORMS FOR 
THE TRAPATT MODE IN A PIN DIODE WITH A DRIVE 
CURRENT USING TEN HARMONICS OF A WAVEFORM 
REPORTED BY SNAPP, et al. 7 


possible explanations for this difference. Firsts we assumed a 10- 
micron-wide depletion layer whereas the depletion width of the experi- 
mental diode was reported to be less. Second, the electrons and holes 
were assumed to have equal ionization coefficients, which should lead 
to the generation of more charge and a corresponding increased period 
of oscillation. Third, we chose a small thermal generation rate whereas 
the experimental diode may have been at a higher temperature. 
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Table 2-9 


NORMALIZED PERIOD T, POWERS , P., AND IMPEDANCES, Z., 

FOR A PIN AVALANCHE DIODE OPERATING IN THE TRAPATT MODE 
= 1.11(0.5 + 0.52 sin (9-6.61°) + 0.29 sin (20-13.9°) 

+ 0.20 sin (36-56.7°) + 0.037 sin (40+207.9°) + 0.041 sin (50+259°) 
+ 0.042 sin (60-86.3°) + 0.028 sin (70 + 2 5 9.4°) + 0.015 sin (80+79°) 
+ 0.008 sin (90 + 239°) + 0.007 sin (100-61°)] 


J = 2 X 10 
s 


Parameter 

Three Harmonics 

Ten Harmonics 

T 

24.31 

26. 56 

P 

o 

0.365 

0.351 

P 

1 

-0.167 

-0.162 

p 

P 2 

-0.018 

-0.023 

P 

3 

+0.00008 

+0.002 

Z 

o 

1.183 

1.128 

Z 

1 

-0.987 + j 0.558 

-0.961 + j 0.582 

Z 2 

-0.343 - j 0.044 

-0.414 - j 0.065 

Z 

0.003 + j 0.788 

0.065 + j 0.795 



The significance of our theoretical results are that we can predict , 
at least qual itatively, the performance of avalanche diodes oscillating 
in the TRAPATT mode for arbitrary current waveforms. 
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3. ANALYTICAL SOLUTION FOR A SECOND HIGH-EFFICIENCY MODE 


A. Description of a Second High-Efficiency Mode 

There are modes of oscillation other than the TRAPATT mode des- 
cribed in Section 2 which provide efficient dc-to-RF energy conversion. 
Figure 3-1 shows a set of voltage and current waveforms that are re- 
resentative of one of these modes of oscillation. In this section an 
approximate analytical solution for the voltage waveform for this mode 
is derived for a rectangular-wave current. First, the essential 
characteristics of this mode of operation are briefly described. 



FIGURE 3-1 VOLTAGE WAVEFORM OF A HIGH-EFFICIENCY MODE OF 
AN AVALANCHE DIODE FOR A RECTANGULAR DRIVE 
CURRENT 
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This mode is distinguished from the TRAPATT mode described in 
Section 2 by the fact that at the beginning of a cycle the diode voltage 
is above avalanche and the external current is assumed to be zero, whereas 
in the ordinary TRAPATT mode the diode voltage is initially high but 
below avalanche and the external current is charging up the diode. At 
Point A in Figure 3-1, the voltage is above breakdown, but there is very 
little particle current in the depletion layer. The charge density builds 
up owing to the avalanche process and drifts out, thereby depressing the 
internal E-field and the terminal voltage. Since the initial internal 
charge is so small, considerable time is required to generate sufficient 
charge to depress the E-field. At Point B enough charge has built up 
that the ionization process acts very rapidly, and a large amount of 
charge is generated during the short time between Points B and C. At 
Point C the field has been depressed sufficiently to stop further 
ionization of charge carriers. We assume that at this point the 
external current begins to flow. The exact instant of time is somewhat 
arbitrary and depends upon the external circuit. For the present dis- 
cussion it is sufficient to assume that the external current comes on 

when the E-field drops to the critical value E below which the avalanche 

c 

process stops. 

As more charge drifts out the ends of the depletion region, the 
field in the center is further depressed. The voltage drops to Point D, 
and the plasma is " trapped” in the middle of the depletion region. This 
portion of the cycle and the portions over to Point F are the same as 
the corresponding portions in the ordinary TRAPATT mode if the plasma 
density is sufficiently large. Thus from Point D to Point E the plasma 
is extracted by the external current, and several transit times will be 
required for large charge densities. 
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At Point E the plasma has been extracted, leaving a residual charge 
of electrons in one half of the depletion layer and a residual charge of 
holes in the other half. As the residual charge is removed, the voltage 
increases more or less linearly from Point E to Point F. At Point F all 
the charge that has been generated internally has been removed from the 
depletion layer. This charge must be greater than or equal to that 
supplied by the external current or the voltage will exceed that at 
Point A. For analysis it is convenient to assume that Point F lies below 
the critical field required for avalanche. Then all the charges except 
those due to thermal generation will have been swept out before the 
diode is charged back up to Point G. The assumption that only a back- 
ground of thermally generated charges exists at Point F enables one to 
calculate the charges generated in going from Point F to Point G. The 
charge in the depletion layer at Point G must be the same as that at 
Point A to have a periodic solution. 

As in the case of our analysis of the ordinary TRAPATT mode, we 
assume that we . have a, PIN diode with equal ionization coefficients, a, 
and equal velocities, v, for both electrons and holes. We showed in 
Section 2 that, when the charge density is small or when the depletion 
region is filled with a plasma, the E-field is effectively uniform. 

These conditions are satisfied over certain portions of a cycle of the 
high-efficiency mode discussed here. 

B . Initial Avalanche and Generation Phase 

During that portion of the cycle from Point A to Point B in Figure 
3-1 the internal charge is small, and there is no external current. The 
operation of the diode over this portion of the cycle is described by 
Eqs . (2-9) and (2-10). Between Points B and C a large plasma is 
generated by the avalanche process. As the charges in this plasma 
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drift, they depress the E-field. However, the avalanche process is 

essentially uniform across the depletion layer, and the E-field remains 

relatively uniform. We assume as a first approximation that E remains 

uniform between Points A and C. At Point B, Q has increased from a 

small value to a magnitude of about 0.1; during this time E is nearly 

constant and equal to E . Similarly, a is equal to a , a maximum 

m m 

value corresponding to E . Between Point B and Point C, Q increases 

m 

rapidly. The E-field begins to decrease, forcing a to decrease. Since 

the mathematics becomes very complicated for a not a constant, we let 

a = a until E reaches a critical value, E , at Point C. For values of 
m c 

E less than E , we assume a = 0. This assumption is similar to the one 
c 

used in analytical solutions for the ordinary TRAPATT mode. The validity 

of these assumptions and the value of E can be established from a computer 

c 

program such as UNSAT. 

Using the above assumptions, we can write the equations for the 
charge density and the electric field during the initial avalanche and 
generation phases of the high-efficiency mode: 


dE 

dt 


-v c 


Q 


(3-1) 


2v (o' - l)t 

Q = Q e s m 
A 


(3-2) 


where Q is the initial charge density at Point A, v is the saturated 
A s 

velocity and o' is the ionization coefficient corresponding to the 
m 

maximum field at Point A. Substituting Eq. (3-2) into Eq. (3-1) and 
integrating, we obtain 


E<t) ■ e a - — to -1) 


2v (a - l)t 
s m 

e - 1 


s m 


0 < t < t , (3-3) 

g 
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where J = v Q is the small particle current at the beginning of the 
A s A 

cycle, t is the time when E drops to the critical field E , and the 
g c 

generation phase is complete. We solve for t from the above equation 

g 

by setting E = E : 
c 


t 

g 


1 

2v (o' - 1) 

s m 


1 + 


2v (a - 1) 
s m 



(3-4) 


The time for the generation phase is a function of both E and J . Since 

A A 

can vary over several orders of magnitude , it can have a dominant 

effect on the generation time. The voltage over this portion of the 

cycle is just equal to E(t), since we have assumed a uniform electric 

field. By substituting the value of t back into Eq. (3-2), we can 

g 

calculate the magnitude of the plasma charge, Q , in the depletion region 

P 

at the end of the generation phase: 


Q = Q A + - 1) (E - E ) 

P A m Ac 


(3-5) 


The plasma charge is a function of only the ionization coefficient 
and the difference between the magnitude of the E-field at Point A and 
that at Point C. 

By analogy with Eqs. (2-20) and (2-21) , the total positive or nega- 
tive charge generated in one cycle, G, is equal to the following: 

G = o' (E - E ) + i Q . (3-6) 

m A c 2 A 


The integral of the external current over one period must equal G, the 
total charge generated internally, for the oscillations to be periodic. 
This sets a limit on the maximum value of J to prevent a second avalanche 
within the same period. 
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We can rewrite Eq. (3-5) for the total charge at Point C as a 
function of the generation: 

Q = 2G - 2 (E - E ) . (3-7) 

p Ac 

As in the ordinary TRAPATT mode, the plasma charge density is equal to 
the total charge generated less those charges that have drifted out of 
the depletion layer. It is the drifting charges that actually depress 
the E-field and stop further ionization. 

C. Field Depression and Charge Removal 

The operation of the device over those portions of a cycle from 

Point C to Point F in this high-efficiency mode is the same as that in 

the ordinary TRAPATT mode described in Section 2 if Q » J. The dis- 

P 

cussions and corresponding equations are the same if Eq . (2-22) in 

Section 2 is replaced by Eq. (3-7) above. If Q < J, then different 

p ~ 

expressions for the voltage during the field-depression and charge- 
removal phases of a cycle must be derived. Since the latter inequality 
corresponds to a less efficient mode, it is not considered in this 
report. 

D. Final Charge-up Phase 

If the drive current has not exceeded the value that would cause 

second avalanche, all the charge will be removed from the depletion 

region at Point F (Figure 3-1). The E-field will again be uniform 

across the sample. From Eq . (2-42) the voltage at Point F is V : 

F 

^F E p + 2 (v + v ) * (3-8) 

s p 
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The diode charges up like a linear capacitor until the field reaches the 

point where significant avalanche occurs, i.e., where the normalized 

ionization coefficient approaches unity. During this time the charge 

density will remain at a small value , determined by the diode back- 

biased saturation current, J : 

s 

J 

Q = — « 1 . (3-S) 

S V 

s 


The external current will continue to drive the field higher into 
avalanche until the voltage reaches the value at Point G. Since the 
charge density is small throughout this portion of the cycle, the time 
required for the final charge-up phase is given by the following: 


t = (E - E )/ J 
c A F 


(3-10) 


where E 

A 


E owing to periodicity. 


To calculate the charge, Q , we must assume a particular functional 
dependence on E for oi • We choose an exponential function, for mathe- 
matical convenience, as done in Section 2. 



(3-11) 


Since we have assumed that the E-field increases linearly with the external 
current, we can solve for the charge from Eq . (2-9) during the final 

charge-up phase of the cycle: 
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(3-12) 


Q(t) = Q 

s 


exp 


2v a 
s 1 

f XJt - ll 

- 2v t 

XJ 

l e J 

s 


where is the value of the ionization coefficient when avalanche 
starts and hence is equal to one; is the charge due to thermal genera- 
tion and corresponds to the diode saturation current. To find Q we 

A 

must substitute for t , the time required for the field to increase from 

E to E • where E is the value corresponding to o' =1. Thus we have 
1 A 1 1 

t = (E - E )/J and 
A 1 


Q = Q exp 
A s 


2v O' 

1 HE. - E,)| 

2v | 

s 1 

XJ 

le A 1 1 

s \ 

- — (E - E) 
J A 1 


(3-13) 


For a self-consistent solution, should be used in Eq . (3-3) for the 

electric field during the generation phase of the cycle. The only 

quantity in our semianalytical solution that is arbitrary is E in 

c 

Eqs. (3-4), (3-5), and (3-6). As mentioned earlier, we can determine 
E^ by comparing our results for the generation as computed by Eq. (3-6) 
with those obtained by using a computer program such as UNSAT. 

No specific calculations of diode parameters for this mode of 
operation have been computed as yet. We can, however, make a few general 
comments. The higher the diode is driven above avalanche before the 
diode current is turned off, the shorter the time before the subsequent 
plasma is created, the field depressed, and the external current 
increased. This results in a rectangular current waveform rather than 
a square wave and decreases the magnitude of the fundamental component 
of current. On the other hand, if the diode is not driven very far into 
avalanche, the density of the generated plasma will be small. A 
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smaller plasma results in a shorter period of oscillation and tends to 
decrease the efficiency. But the duration of the high-voltage state 
is increased, resulting in a more symmetrical waveform, which tends to 
increase the fundamental efficiency. The interaction of these various 
trade-offs is not clearly understood at this time. Additional investi- 
gation of this mode of oscillation is required. This mode may be more 
representative of the high-efficiency oscillations observed experimentally 
where the fundamental frequency is 1/2, 1/3, or 1/4 the transit time 
frequency. 


53 



PRECEDING PAGE BLANK NOT FILMED 


4. TRAPATT OSCILLATOR WITH MULTIPLE AVALANCHE PER CYCLE 
A. Introduction 

In this section the behavior of a TRAPATT PIN diode is analyzed 
under the condition of multiple avalanche per cycle of oscillation. 
Previous analysis of TRAPATT oscillations have avoided this possibility 
by assuming the diode current to be below a certain threshold value 
(see Section 2 and Refs. 3 and 4). The results of this section are an 
approximate solution for the voltage and current waveforms over an entire 
cycle of oscillation. The problem considered can be summarized as 
follows. A dense plasma is formed uniformly throughout the depletion 
region of the PIN diode by avalanche ionization, as in the ordinary 
TRAPATT mode. Part of the trapped plasma is removed near the ends of 
the depletion layer by the external current, leaving a residual charge 
of electrons at the left side of the plasma and a residual hole con- 
centration on the right side. The residual charge outside the plasma 
causes the electric field to increase toward either boundary. Since 
the magnitude of the residual charge, and hence the boundary value of 
the electric field, is proportional to the external current during this 
segment of the cycle, the E-field near the edges will exceed the value 
for avalanche breakdown for sufficiently large currents [see Eq. (2-38)]. 
The electrons created at the left boundary by this second avalanche 
exit, but the holes drift to the right toward the slowly receding plasma, 
and increase the residual charge density in that region. Similarly, 
holes created at the right boundary by the second avalanche exit, and the 
electrons drift to the left toward the plasma. 
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The dynamics of this second avalanche, and any subsequent avalanche, 
are not completely understood at this time. We will assume that the 

transients associated with this avalanche are short lived (on the order 

/ 

of one transit time) and allow the diode to reach a new quasi-steady 
state. The new field configuration and charge distributions are 
determined from a quasi-static approximation to the continuity equations 
and an assumed constant external current that exists over the majority 
of the first half of the cycle. This current decreases near mid-cycle 
to hold the voltage just below uniform avalanche and then drops to zero 
over the last half cycle. The quasi-static approximation is based upon 
the assumption that the only time variations in electron and hole dis- 
tributions are caused by the slowly moving boundaries of the trapped 
plasma. Since the plasma velocity is much less than the saturated 
velocity, this approximation seems reasonable, 

B. Initial Charge-up and Plasma Generation 

We shall assume that initially a PIN diode is biased just below 
avalanche with a uniform electric field, E^, over the depletion layer 
and a constant current, J, applied to the device. The form of the 
solution is the same as given before in Eq. (2-11); 

V(t) — V + J t , 0 t 5 t . (4-1) 

A m 

As the current increases, the rate of voltage growth increases, but the 

time over which this formula applies, t , decreases. The peak value 

m 

that the voltage attains increases about two fold for a ten-fold 
increase in drive current. 

The plasma generation sequence then follows in which the carriers 
generated by the avalanche reduce the field across the depletion region 
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rapidly from its maximum value E 

m 

plasma. The total charge filling 
Eq. (2-19) 


to the value E within the dense 
P 

the depletion layer Q is given by 

P 


Q = Q + - 1) [E - E + Jt ] , (4-2) 

pm m meg 


where Q = J/v , v = the saturated velocity, a is the ionization 
mss m 

coefficient corresponding to the field E , and t is the time of this 

m g 

process given by Eq . (2-18). The voltage drops rapidly from the maximum 


E to E as the dense plasma removal begins. The velocity at which the 
m p 

dense plasma is removed is the plasma velocity v = J/Q correspondong 

P P 

to the low plasma field, E . So far the diode behavior is similar to 

P 

the lower current case in the ordinary TRAPATT mode. 


As the dense plasma is removed, its boundary moves toward the 

<? 

center of the depletion region at v , and the electric field increases 

P 

in the region outside the plasma. This continues until the field at 


the boundary reaches the avalanche value E . In the region of residual 

c 

charge the field is given by 


E(x) = E + n (x -x) , 0 < x ^ x , (4-3) 

p r p p 


where x = v t is the position of the plasma boundary and n is the 
P P r 

residual charge density. The time, t , it takes for the field at the 

c 

edge to reach the avalanche value is given as 


E - E 

t = — (v + v ) . (4-4) 

c v J s p 

P 


At this time the second avalanche is established near the depletion 


layer's edge. To determine the field when this occurs, we invoke a 
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quasi-static approximation. The normalized static continuity equations 
for static fields can be combined with Poisson’s equation (neglecting 
diffusion) to yield a second-order nonlinear equation for E: 


v 



2 cv J 



(4-5) 


Here we have also assumed equal hole, electron ionization rates, ot y and 
equal carrier velocities, v. 

J, the external current, equals J + J , the sum of electron and 

n P 

hole currents, and v 7 equals ©V/3E). The boundary conditions on the 
left side of the junction are 


9E 


ax 


x=0 


J 


(4-6) 


and 


E = E at x - x 
P P 


(4-7) 


The field profile will be symmetrical about the middle of the depletion 
region, with the E-field continuous throughout this zone. 

Since both a and v are functions of E, one method of solution for 

Eqs. (4-5), (4-6), and (4-7) is to separate the problem into two regions. 

In one region the field is above the critical field E for avalanche 

c 

and the other below, so cy may be neglected. Near the boundary, in the 

second avalanche zone v equals v and Eq. (4-5) reduces to 

s 


2 

a e 
2 

a x 


= 2 a J 


(4-8) 
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over the range 0 ^ x ^ x , where x is the position at which the field 

c c 

equals E . The power representation of the ionization coefficient 
c 

dependence on field will be assumed 



(4-9) 


as suggested by McKay, where oi is a constants 0.2 for silicon. This 

o 

form of the ionization coefficient is used rather than the exponential 
for mathematical convenience. Equation (4-8) can be written as 



2 & J 
o 

6 

E 

c 


(4-10) 


and has a solution given by 


E(x) 


1 

1/2 

(2Ax + 2C) 


0 < x £ 


where 



> 


(4-11) 


and using the boundary condition 



c 
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Setting E(x) = E and solving for the position x yields 
c c 


x 

c 



(4-12) 


The region over which second avalanche occurs (0 £ x < x ) is not a 

c 

strong function of x being about 0.22 for currents just large enough 

c 

to cause second avalanche; x decreases to about 0.18 at four times 

c 

this threshold current. 


The electric field at the boundary is 


E (0) = E 

c 



1/6 


(4-13) 


and also is seen to be a slowly varying function of the current. The 
voltage across the second avalanche region is quickly established, in 
a fraction of a transit time, and then remains constant so long as the 
extracted current is constant. The time-dependent voltage is 


V(t) 


(2Ax + C) 1/2 

c 

A 



t 

t 

a 


t 

V — , 0 £ t £ t , 
s t a 

a 


(4-14) 


where t is the time for a carrier to traverse the second avalanche zone 
a 

= x c / v s j and A and C were given above. Between the second avalanche 

region (ending at x = x ) and the dense plasma region (starting at 

x = x ), there is a region sufficiently below avalanche that a in 
P 

Eq. (4-8) is neglectable, but of high enough field that v^ v over 

s 

most of its range. Equation (4-5) thus reduces to 
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0 


(4-15) 


whose solution is linear in space 


E(x) = E -I (E-E) , x ^ x ^ x . (4-16) 

c x - x cp c p 

L p c J 


The voltage across this region is 


/e - E \ /x 2 x 2 

V(t) = Ex - Ex -I- 2 I j-^ - -£■ 

cp p c 1 x^ - x c / I 2 2 


(4-17) 


Since x = v t. x = constant, and E » E , we have 
p p c c p 


V(t) as V + — (v t - x ) , t fits . 

s 2 p c c v 

P 


(4-18) 


Thus the total voltage for the plasma extraction segment under second 
avalanche is the sum of the voltage across each region, given by 
Eqs . (4-14) and (4-18): 


V(t) =V + — (v t - x ) , 

s 2 p c 


(4-19) 


where 


c 0.5 

t + — £ t £ 

c /v v 

s p 


and is linear with time over this range. 
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The final step is to remove the sparse plasma left behind from the 

original avalanche, and from the continuing second avalanche. The field 

within the second avalanche region is a complex quantity somewhat 

different than in the region x ^ x ^ 0.5. Analysis has shown, however, 

c 

that a linear approximation to the field E(x) over the entire depletion 

region (this approximation accounts for the extra charge generated within 

the second avalanche and matches the slope at x = x ) yields reasonable 

c 

results for the excess charge removal phase. Figure 4-1 shows the 

normalized field distribution in the presence of second avalanche. The 

* 

total excess charge is n ^ (E(0) - E )/0.5 and the field throughout 

r p 

the depletion region is 

E (x) = E + n (0.5 - x) , 0 £ x ^ x , (4-20) 

p r r 


where x is the extent of the residual charge measured from the left 
r 

boundary. The value of x decreases linearly with time: 

x = 0.5 - vt , 0 <. t ^ t , (4-21) 

r s i 


where t will be given below. The corresponding terminal voltage is 

E (0) - E 

V(t) = E + (1 + 2v t) , 0 £ t £ t . (4-22) 

p 2 s i 


This phase of the cycle continues until the terminal voltage reaches 

its initial value. The current must be decreased at this time to prevent 

the voltage from overshooting its initial value. Hence we will require 

that J decrease in such a way that V(t) is equal to V over the final 

A 

segment of the cycle. The time t corresponds to the voltage reaching 

its initial value V , 

A 

62 




63 



(4-23) 


2E - E - E 
t = c ° P 

i 2E v 

o s 


At this time the area below E(x) and E to the left of x must equal 

c c 

the area above the E(x) and E to the right, until x = 0. 5. From the 

c 

time t until t = 0.5/v the current decreases at a rate that keeps 
is 

the two areas equal (see Figure 4-1), and the current time dependence 
will be 


J(t) = J 




2 Qi 

o 

J 


0. 5 
v 

s 



t £ t £ 
i 


0.5 

v 

s 


(4-24) 


It is seen that this portion of the cycle lasts only a fraction of a 
transit time, and at its conclusion the current is reduced to zero, which 
maintains the voltage at its high initial value ready for the start of 
the next cycle. 

Multiple avalanche thus has two effects on the total device per- 
formance. The power output increases almost linearly with increasing 
drive current whereas the dc-to-RF efficiency decreases only slightly. 

For example, calculations indicate a relative decrease in efficiency of 
10 percent for a four-fold increase in current, over the threshold value 
for second avalanche, and a corresponding four-fold power output increase. 
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5. ANALYSIS OF THE UNIFORM, SYMMETRIC, AVALANCHE DIODE 
IN THE HIGH-FIELD STEP -I ONI Z ATI ON APPROXIMATION 


A. Introduction 

A qualitative understanding of the transient and oscillatory 
behavior of a PIN diode as it moves in and out of avalanche ionization 
is sought by approximating the variation of ionization with field as a 
step from none to a fixed level of ionization at a critical field value. 
This critical value is assumed high enough so that electrons and holes 
move at saturated velocities. 

In this section it is shown that a solution can be computed when a 
pattern of regions with field values above and below critical has been 
prescribed. A rule governing admissibility of patterns has been derived. 
The rule has shown that many suggested patterns are inadmissible; one 
family of useful admissible patterns has been found. 

In the approximation, electrons and holes move at the same constant 
speed, electrons to the left and holes to the right, according to assumed 
sign conventions. In the (x,t) plane the electron and particle paths are 
straight characteristic lines. The density of electrons, n, has the 
small thermal value e at the right edge of the diode, is continuous along 
a characteristic running to the left, is constant in regions with no 
ionization, and increases in regions with ionization. Hole density, p, 
satisfies the same conditions on each right-running characteristic. 

The requirement that the field be above its critical value in 
regions with ionization, and below in regions without, leads to the 
following rule: 
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If a boundary segment between regions in the (x,t) plane with 
and without ionization is horizontal, a line of constant t, 
then p = n along the segment. If the boundary segment is 
not horizontal and the region with ionization is to the 
left, then near the segment in both regions p < n. If the 
region with ionization is to the right, p > n. 

The approximate analysis is used to indicate possible patterns of 
avalanche regions in the space-time plane for the simplest case of a 
uniform PIN diode. Qualitative results for the variation of current 
and voltage drop across the diode are deduced. 


B. Basic Equations 

The basic equations are those presented in Section2 in the nor- 
malized form they read 


n 

t 


+ (v n) 
n x 


+ (v p) 

P X 


E 

x 


E 

t 



=J+vn-vp, 
n p 


(5-1) 


where E is the field strength, n is electron density, p is hole density, 

v n electron velocity, v is hole velocity, a and R are ionization 

P 

rates, and are donor and acceptor densities, and J is the external 
current density. Subscripts x and t denote partial differentiation. 

The boundary conditions become 
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p = e at x = - | , 

. (5-2) 

n = e at x = £ , 

where e is the small charge due to thermal generation and the origin has 
been taken for convenience in the central plane of the diode. Note that 
in Sections 2 and 3 the origin was chosen at the left-hand edge of the 
depletion region. 

Introduce function F defined by 

x t 

F = E - P (N - N ) dx - P J dt , (5-3) 

•J x D A J t 

0 0 


where t is a convenient fixed time and x„ is a fixed location. One 
0 0 

finds 


F = p - n , f = v n - v p 

x t n p 


so that 

n = (vF + F ) / (v - v) , 
p x t n p 

(5-4) 

p = (v F + F ) / (v -v) 
n x t n p 


The system is represented by the single differential equation 




v F + F 

p x t 

V - V 

n P 


+ 3 



v F + F 
n x t 

v - v 
n p 

(5-5) 
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with boundary conditions 


v F + F 

P X t T 

= 0 at x = - 5 - , 

v - v * 

n p 

(5-6) 


v F 4 - F 

n x t t 

= 0 at x = - ? 

v - v 
n p 


Since only derivatives of F appear in the differential equation and 
boundary conditions, if F is a solution so is F + C a solution for any 
constant, C. 


C. High-Field Step-Ionization Approximation 


To obtain qualitative understanding of the avalanche process, we 
have simplified the equations as follows. Near avalanche the field will 


be so high that electron and hole velocities are essentially saturated. 
In the normalized units, this means, 


v = - v = 1 . (5-7) 

P n 


The ionization rates are taken to be equal. Their sharp variation from 

a low value to a high one at a critical value of E is taken to be a step 

change from zero to a constant value a . The step is taken to occur at 

E = 1 ; this corresponds to the critical value E used in Sections 2 and 

c 

3 : 


a = 3 


0 


Otr 


E < 1 
E > 1 


(5-8) 
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The resulting approximate equations for F are the wave equation 


when E < 1 and a variant of the telegrapher’s equation when E > 1 


F - F = 0 
tt xx 


for E < 1 , 


F - F =2 cv F for E > 1 

tt xx 0 t 


with boundary conditions 


(5-9) 


One finds 



F + F = 

t X 

-*2e at x = i , 


F - F = 

“2s at x = - l 


t X 


characteristic coordinates are 


r = t + x 

) s = t - X 

F 

= F - F 

ii 

+ 

X 

r s 

t r 

F 

= (Q/2)(F 

+ F ) 

rs 

r 

s 

F 

= - 6 

at r - s = 1 

r 



F 

= - e 

at r - s = - : 


(5-10) 


(5-11) 


(5-12) 


Since 


n=- F r , 


P = - F s ' 


(5-13) 
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useful solutions of F have negative partial derivatives. From the 

general theory for hyperbolic equations, one knows that F is continuous 

with continuous first partial derivatives, except that F and n may have 

r 

a jump in value across a characteristic line s = constant and F and p 


may jump at an r = constant line. 


s 


D . Solution Forms 

Where E < 1, function F satisfies the wave equation and therefore 
is given by 


F = f (r) + g ( s) (5-14) 

for arbitrary functions f and g. These functions are to be chosen so 
that F is continuous and fits the boundary conditions. 

Where E > 1, function F satisfies the telegrapher's equation. The 
value of F at a point :: (r^s^) is ^ ouncl ky appropriate integrations 

involving the Riemann function R, which depends on P and a variable 
point (r,s). In the present case R is only a function of the two 
variables (r - r Q ) and (s - s q ) : 

R = R(r - r Q , s - s Q ) (5-15) 


with 


- o Q (y + z)/2 

R(y,z) = e I Q (<* 0 /~yI) 

where I is the modified Bessel function of the first kind. 

In particular, if F is known along two intersecting characteristic 
lines, r = constant and s = constant, as AC and BC in Figure 5-1, F 
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FIGURE 5-1 CHARACTERISTIC INITIAL-VALUE PROBLEM 


can be found at P from the formula 


F(P) = F(A)R(A) + F(B )R(B ) - F(C)R(C) 

(5-16) 

A B 

- P F(R + a R/2) ds - P F(R + a R/2) dr 
z 0 J y 0 

C C 


Similar formulas hold when the data are given in other ways. 

Two solutions are F > and F : 

1 2 


<* o (r + s) 2 a t 
F = e = e 


F 2 =e 


a Q (r+s)/2 


<y t 
0 


l 0 (a 0 ^Fs)= e ^ xL 7t 2 - x 2 y 


(5-17) 


(5-18) 


In a useful solution, with neither n nor p negative, 
appear with a negative coefficient. 


F or F would 
1 2 


E. Avalanche Regions 

A solution to the equations may be computed when initial conditions, 
Uf D ~ n a ^ and J are given. For simplicity, the case of a PIN diode, 
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with equal donors and acceptors, = 0, is considered here. By 

Eq. (5-3) we have, 

t 

E = F + Q , Qj= J J dt . (5-19) 

t 

0 

One way of specifying J or Q is to assign regions of avalanche and 
no avalanche in the (x,t) plane. For then F can be computed by solving 
the differential equation and boundary conditions in set of Eqs. (5-12). 
And, since E = 1 on the boundaries separating the regions, Q can be 
found; and this gives E everywhere. 

Let a region with avalanche ionization be an A-region and one with- 
out, a W-region, the letter W being a reminder that the wave equation is 
controlling in such regions. Since holes and electrons are neither 
created nor destroyed in a W-region, the value of p is constant along 
each characteristic running to the right; the value of n is constant 
along a characteristic running to the left. In an A-region the value of 
p is increased by the ionization as one follows a right-running charac- 
teristic, while n increases along those running to the left. 

Consider a segment of boundary between an A-region and a W-region. 

On the segment, E equals 1. If the segment has slope zero in the (x,t) 
plane--i.e., follows a line of constant t — then Q and F are constant 
along it. If the segment has nonzero slope, it is crossed by lines of 
constant t. On such a traversal in the A-region E falls as it approaches 
the boundary and continues to decrease as one leaves the boundary in the 
W-region. Thus, whether there are discontinuities in derivatives at the 
boundary or not, E has a negative sign on both sides of a boundary with 
the A-region to the left; it has a positive sign if the A-region is to 
the right. By Poisson’s equation, the third equation in set (5-1), and 
Eq. (5-13), the following rule is established: 
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n < n or F < F , if the A-region is to the left, 

F r s 

p > n or F < F , if the A-region is to the right, 

s r 

p - n or F = F , if the boundary has constant t. 

r s 

The condition at the left edge of the diode, x = - -g* , is p = e. If 
there is a W-region adjacent to this edge, then p = e along any charac- 
teristic starting from that edge. Where the characteristic reaches a 
boundary with an A-region, in which n will be greater than e , by the rule 
above, the A-region is to be to the left. Since the characteristic line 
has slope 1, the boundary has to have zero or positive slope less than 
one. The same restriction holds with a change of sign at the right 
edge. In more graphic terms, as time increases, the boundary of a W- 
region adjacent to an edge of the depletion layer overtakes character- 
istics issuing from that edge in the W-region. 

F. The Characteristic Pattern 

If F and F and their first derivatives with x are continuous at 
t 

t = 0 across the diode, - ^ ^ x ^ , then any discontinuities in solution 

can only be along characteristic lines originating at the edges of the 
diode together with their reflections from the edges. These character- 
istics are shown in the sketch. Figure 5-2. It seems reasonable to 
seek a solution in which these characteristic lines separate the A- 
regions from the W-regions. However, this attempt fails, as the following 
arguments show. 

If triangle A C A^ is a W-region, the whole strip A A B B is also, 
112 ' 1 2 3 2 9 

since none of the assumed boundaries are overtaking. For the same reason, 

strip B B A A is also a W-region. If triangle AC A is a W-region, the 
^ u 4 O Z Z o 
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FIGURE 5-2 POSSIBLE LINES OF DISCONTINUITY IN AN INITIALLY 
SMOOTH SOLUTION 


whole diode remains without ionization, a case of no interest here. If 
triangle AC A is an A-region, it is to the left of a W-region in which 

/ ^ O 

n = g • This violates the rule. 


Since the various triangles are symmetrically placed, the argument 

of the above paragraph shows that no triangle can be a W-region. The 

possibility that quadrilateral C B C A is a W-region is now explored. 

12 2 2 


Triangle A C A is an A-region in which the hole density p is e 
112 

along edge A A , and increases along boundary AC to a larger positive 
12 11 

value at . The hole density increases in the same way along other 
characteristics parallel to A^C_^ in the triangle but for shorter periods 
of time. Thus along C A . p has a positive value at C and decreases to 
e at A . These values are carried across the quadrilateral W-region 


74 


unchanged so that p decreases from a positive value at B g to e at c 2 * 
Along the characteristic in the triangular A-region B 2 C 2 B 3 > n 

increases from e at B^ to a larger positive value at C^. Thus, near C g , 
n is positive and reasonably large, whereas p is near e. This violates 
the rule p > n when the A-region B 2 C 2 B 3 is to the ri g ht - Thus no solu- 
tions exist with only the characteristic lines of the figure as region 
boundaries. 

G. An Admissible Pattern 

A pattern satisfying the rules is shown in Figure 5-3. Time inter- 
vals a and b are variable on their ranges: 0 ^ a £ O^b. The whole 

pattern is periodic with period (2 + b) transit times. The type of each 
region is marked on the figure. 
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FIGURE 5-3 ADMISSIBLE PATTERN 
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On constant time line B C , F is constant and p = n = e . One may 

o o 

readily check that the left-running characteristics to points near 

boundary D G M P have a longer run in A-regions than the right-running 
o o o o 

ones do. Thus, n is greater than p along this boundary with an A-region 
to the left in accordance with the rule. The rule is also observed on 


D H N Q by symmetry. In the remainder of the time, the electrons and 
O O o o 

holes are swept away so that their densities are again near zero on 


B C . 


1 1 


The pattern can repeat. 


The two characteristics reaching a point in triangle BCD have 

O o o 

equal length back to their crossings of B C , where p or n is zero. Thus 

o o 

p equals n in the triangle and E is zero. On each constant t line 

x 

cutting the traingle, E is constant across the triangle and increases 
outside the triangle at both ends toward the diode edges, as in Figure 
5-4. The value of E on x = 0, E(0,t) is determined from E(0,)) = 1 and 
by the last of Eqs. (5-1) 


= J - (p + n) . (5-20) 

There is great freedom in choosing J so that E(0,0) = 1 and E(0 ,tt) =1, 
and E (0, t) > 1, for 0 < t <■§■. 

x 

T A-8327-34 

FIGURE 5-4 FIELD VARIATION ACROSS THE DIODE, 0 < t < y 2 
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Current J is fixed t^l + bby the requirement that E = 1 

on the boundaries between A-regions and W-regions. At a typical value 
of t on^ t ^ 1 + b - a, the variation of E would be as in Figure 5-5. 
On 1 + b - a ^ t ^ 1 + b, the points at which E = 1 move to x = ± b • 



FIGURE 5-5 FIELD VARIATION ACROSS THE DIODE, % < t < 1 + b - a 

In the interval , l + b^t^2+b, maximum values of E occur at 

each time at the diode edges with some decrease in value to a minimum 

in the center. During the last half transit time, E is constant at 

each time across the central section in triangle R C B . Here also J 

oil 

may be chosen quite arbitrarily except that E < 1 in the time interval, 
and E = 1 at t = 2 + b. 

H. Analytic Solution of a Special Case 

At time zero, only thermally excited electrons and holes are 
present. The initial condition is 


p - n = e 


with e a small constant. For function F, we have, 


F = 0 , F = - 2e , t = 0 , 
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where the initial constant value of F has been taken to be zero. 

The solution in the initial A-region, B C L D K B of Figure 5-3, 

o o o o o o 

is the same as it would be if the diode were in avalanche for all time. 
This suggests applying Laplace transform techniques to obtain the 
solution. 

The solution process is simplified if function F is replaced by G: 

-at 

G = e ° F . (5-21) 


The differential equation and side conditions for G read: 


2 

G - G = a G 
tt XX o 




G = 0 , G = - 2 e 


at t = 0 , 


- a t 

G + g + a G = -2ee ° at t = k 

t x o 


(5-22) 


G = 0 at x = 0 

x 


The Laplace transform of G is G: 


G = J e ^ G dt 
0 


(5-23) 


Conditions on G become 
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2 - 2 2 2 
G — q G - 2 e , q — |JL ~~ Of 
xx . o 


nG + G + Ct G 
x o 


u + a 


at x = 


G = 0 
x 


at x = 0 


From these one readily finds 


5--H - 

q v 


2a cosh qx 
o 


[<U + a ) cosh (q/2) + q sinh (q/2)] \ 


The first term, -2g q , has the inverse transform 


/ - a t a t\ 
(e/<* o > (e ° - e ° ) , 


so the corresponding part of F is : 


2a t 

F = (e/a ) (1 - e ) 
I o 


Thus, Fj is a multiple of the solution F = F^ found in a previous section 
plus a constant. 

The second term of G may be expanded in a series of exponentials: 


G = jU (G k+ + \- y ’ 


_ 4 ««o / ^V! \ e ~q (k+l+2x)/2 

Gk± = q 2 (p+a +q ) \ ^o + V ° 
o 


79 



It can be shown that the inverse transform of G is zero for times t 

k± 

less than (k + 1 + 2x)/2. Thus, no terms of this series contribute to 


the solution in triangle of B C D of Figure 5-3: F = F in this 

o o o I 

triangle . 


Only the inverse of contributes in C D L , while only the inverse 

0+ o o o 

of G is used in B D K . These inverses have the same form except for 
0- o o o 

the sign of x. Now 



4ea 

o 


q (jjL+a +q) 
o 


-q (l-2x)/2 
e 


may be written 



2g_ 

2 

q 


a 

o , . 

— - (n/q - i) 


-q(l-2x)/2 

e 


(5-28) 


where the first factor and each term in the brackets with the exponential 
function appears in transform tables. The inverse transform is 


G = 2g f sinh qi (t - t / ) C(t') dt' , 
°* 0 


where 


C(t) = H 


i (a y, y 0 ) I, (a y, y_) 

0012 y 1 o 1 2 
z 


y i = Jt - (1 - 2x)/2 , 


(5-29) 


(5-30) 


y 2 =V t + (1 - 2x)/2 , 
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and H is the Heaviside function, which is zero for t less than 
(1 - 2x)/2 , and one for t greater than that value. 


Finally, in triangle C D L , we have. 


F = F + G . (5-31) 

I 0+ 


In triangle B D K , function F has the same form with x replaced by (-x) . 
o o o 


The simplest special case with a = b = 0 is considered here; it is 

sketched in Figure 5-6. On D L , t equals x + %, so s equals ^ . It 

o o 

follows that 


t = (r + p- )/2 , x = (r - %)/2 on D L . (5-32) 

o o 


Let function F, as given in Eq. (5-31) but with these values substituted 
for x and t, be f(r). 



FIGURE 5-6 PATTERN WITH a = b = 0 
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On D K , t + x equals i> so r equals ^ . It follows that 
o o 


t = (i + s)/2 , x (jg* — x)/2 on D K . (5-33) 

o o 


Now function F is given, in triangle B D K , by Eq. (5-31) with x 

o o o 

replaced by (-x) . Hence, on D K , F equals f(s) with the same function 

o o 

f as before. 

In the characteristic quadrilateral D L R K , function F is given 

o o o o 

by 


F = f (r) + f (s) - fft) in D L R K (5-34) 

o o o o 


since r = s = ^ at D . 

o 


The solution in triangle R L C is F = f(s) + g(r), for some 

o o 1 

function g. The boundary condition along L C is F = - e . Hence, 

o 1 r 

g is - er plus a constant chosen so that F will be continuous across 

R L , where r = 3/2 : 
o o 


F = f(s) + f (3/2) - f(i) - e (r - 3/2) in R L C . (5-35) 

o o 1 


By the symmetry, we have 


F = f (r) + f (3/2) - f(i) - e (s - 3/2) in R K B . (5-36) 

o o 1 


Finally, in triangle R B C the solution matching the known values 

oil 

on R B and R C is 
o 1 o 1 


F = 2 f ( 3/ 2 ) - f(|) - e (r + s - 3) in R B C . (5-37) 

oil 
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Current and Voltage 


I . 

The total external current density J is to be chosen so that E > 1 
in A-regions and E < 1 in W-regions. The relation between E and F, 

Eq. (5-19) may be written 

E = F + Q + 1 , Q=[’ t Jdt . (5-38) 

J J ' 0 

This is derived from Eq, (5-19) by noting that parameter t should be 
chosen so that E = 1 and F = 0 at t = 0. 

For t <b, the requirement is that Q + F > 0. Since Q is constant 
45 7 J J 

for each value of t while F varies, this inequality is satisfied by 
taking 


Q > - min F , 0 < t . (5-39) 

J x 

Since E - 1 on K D L , we have, 
o o o 

Q = - f (r) , ^ tsi . (5-40) 

J 

A further requirement is that (F + Q ) is negative inside the angle * 

J 

K D L and positive outside it. The requirement is shown to be satisfied 
o o o 

in two particular cases: Put F = F(x,t). Then one instance of the 

requirement is 

F(0,fc)<F$,fe) 

Substitution from Eq. (5-31) and some simplification give 
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a / 2 

0 < e ° - I (a /2) 

o o 


But this is a known inequality, so the requirement is satisfied. 

If a <4, the whole line t =i satisfies the condition because the 
o 

inequality 

f(o,|) < F(x,b) , 0 < X <; | , 


readily reduces to 
a x 

0<e° - I (a y) + [c] , y = x (1 - x) , 

o o 

where [C] stands for the integral term containing C. This inequality 

holds provided it does so when positive [c] is dropped and I (<y y) is 

* 22 o o 

replaced by the larger quantity, exp (a y /4) : 

o 

2 

ax a x(l-x)/4 
0 < e ° - e ° 

This inequality is satisfied if 


a (1 - x)/4 < 1 
o 

Since 0 < x this requires a <4. Note that a < 4 is sufficient 

o o 

for the requirement to hold, but may not be necessary. 

In the second half of the period, the whole diode is out of 
avalanche. The restriction on the current reads 


* 

Reference 11, Section 2.11. 
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Q < - max F 
J x 


1 < t < 2 


(5-41) 


The voltage drop, V, is the integral of the field. From Eq. (5-20) 
is given by: 

i 

s 

v = j - [ (p + n) dx . (5-42) 

t ‘ “ §■ 

Since p + n = - F^_, this may be integrated over time to obtain 

i 

V - V = Q + M , M=f 1 Fdx , (5-43) 

O J - 2 

where V is the diode voltage at time zero, 
o 

From the analytic solution forms given above, one may compute F and 
M. As a guide to the computations and to an understanding of the solution, 
a qualitative sketch of (-M) is given in Figure 5-7. In Figure 5-7, 
curve OABC is a sketch of the other functions needed on the various time 
intervals: (- min F), (-f), and (- max F). 

The curve for Q in Figure 5-7 would be a continuous curve lying 
J 

above OA, running along AB, and falling below BC. Two extreme cases 

are shown: the dashed curve OABC itself, and the dotted curve, ODABGC , 

departing a "large" distance from OABC. Corresponding sketches of 

(V - V ) and J are given in Figure 5-8. 
o 

If Q follows its boundary curve OABC, the resulting variations 
J 

of (V - V ) and J appear from Figure 5-8 to be roughly in phase, a 
o 

condition of low efficiency. Higher efficiencies are obtainable by 
following an appropriate curve ODABGC with large departure from the 
boundary curve. 
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FIGURE 5-7 AUXILIARY FUNCTIONS 


J . Conclusions 

The approximate analysis of this section has been used to indicate 
possible patterns of avalanche regions in the space-time plane for the 
simplest case, the PIN diode. Estimates for the variation of current and 
voltage drop across the diode have been deduced without computation. 

The limitations on the current have been found, and first steps have been 
taken in selecting a current variation for efficient operation. Further 
efforts to refine the methods of pattern and current selection for the 
PIN diode and to extend them to diodes of other types are warranted. 
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6. EXPERIMENTAL RESULTS 

A . Summary of Experimental Work to Date 

The principal objective of our experiments was to observe the 
effects on efficiency and output power of the various and independently 
variable vector impedances "seen" by an avalanche diode at each of 
several harmonically related frequencies. 

The use of channel -dropping diplexers is appropriate for controlling 
loads that are separate and independent with respect to frequency. 
Accordingly, diplexers that dropped channels centered at 1500, 3000, and 
6000 MHz were constructed and used with an avalanche diode having a break- 
down voltage, V , of approximately 60 volts and an IMPATT frequency, 

B 

f , of about 20 GIz. Experiments with this circuit were discon- 

IMPATT 

tinued after preliminary results yielded only low efficiency at 1500 MHz 
and low output powers. In retrospect, it is seen that the pulsed bias 
supply then in use did not have the proper source impedance and was 
power limited. This experiment should be resumed with an appropriate 
bias supply which is now available. 

A second experimental circuit was arranged for use with inexpensive 

diodes having V ^ 160 volts and known to oscillate very efficiently at 
B 

a fundamental frequency near 750 MHz. The diplexers dropped channels 
centered at 750 and 2250 MHz while the RF circuit also included a stub 
that presented an open circuit to the diode at 1500 and 3000 MHz. The 
system was tried first with a bias supply now known to be inappropriate, 
and then with an almost ideal bias supply. In either case, the "high- 
conductivity state" was approached only feebly, and output efficiencies 
were too low to be of interest. Measurements showed that the circuits 


89 



contained too much loss, which prevented us from obtaining the required 
impedance transformations. The real part of the impedance presented to 
the diode in the desired frequency range was always greater than 10 ohms, 
regardless of tuning. For these diodes the real part of the circuit 
impedance at the fundamental frequency should apparently be less than 
10 ohms when the real part of the circuit impedance at the second har- 
monic is relatively low. 

Near the end of the contract, it was decided to set the diplexers 
temporarily aside and obtain information on the effects of harmonic 
tuning by indirect means. Two types of circuits were constructed, and 
the high-conductivity state (HCS) and accompanying high-efficiency 
oscillations (HEO) were obtained. The construction of these circuits 
allowed the vector impedance of the circuits to be measured and referred 
to the diode junction readily and with high accuracy. A great many 
different combinations of the settings of the tuning elements were used. 
All observations of HCS and HEO corresponded to low values to the real 
part of the circuit impedance for the first few harmonics. All cor- 
responding values of the imaginary part of the circuit impedance were 
observed to be inductive for the first and second harmonics. 

A few observations were also made with series and parallel combina- 
tions of avalanche diodes and with a novel and potentially useful 
arrangement in which a pair of diodes is in series for biasing purposes 
but in parallel with respect to the RF circuit. 

B . Terminology and General Observations 

This subsection lists some terminology for the purpose of simpli- 
fying and clarifying the discussion of the experimental results. Some 
experimental observations are also included that have not previously 
been published and that suggest future areas of investigation. 
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The reverse voltage required to draw a current from 1 to 100 \xA 

was noted as the breakdown voltage, V . In general, we found that 

R 

Fairchild FD-300 diodes that had 160 < V < 180 volts gave the high- 

B 

efficiency oscillations in the frequency range 650-900 MHz. In a batch 

of these diodes, this range of V is the lowest to be found since the 

manufacturer provides that V be not less than 150 volts. When diodes 

B 

were placed in parallel, they were chosen to have, as nearly as possible, 

the same V . 

B 

Reverse bias was applied from a pulse supply that can be considered 
to be an ideal voltage source, V , with a series resistance R . It 

0 LLi 

was found advantageous to have, insofar as possible, R variable. When 

LL 

V > V , a "static breakdown" current, I,, flows and V - V = R I , 
0B 1 0 1 LL 1 

where V is the voltage across the diode junction. The point* (I^V^), 
in Figure 6-1 will be called the "low-conductivity state" (LCS). 



FIGURE 6-1 VOLTAGE AND CURRENT RELATIONS IN AVALANCHE- 
DIODE EXPERIMENTS 
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Certain critical tuning settings of the RF circuit created a 
condition whereby the diode spontaneously and rapidly (typically in 
20 ns) switched to a high-conductivity state (HCS) defined by the 
point, (I , V ) , in Figure 6-1. For FD-300 diodes V was typically 

2* 2t 2 

about 20 percent less than V . The ratio of I to I was typically 

1 2 1 

5 to 1, but ratios from 1.2:1 to 12 :i have been observed. RF oscil- 
lations at low frequencies were obtained for both low-conductivity and 
high-conductivity states. The LCS oscillations were characterized by 
considerable noisiness in the RF pulse envelope whereas the HCS oscil- 
lations as seen in the pulse envelope were relatively free of noise. 

These two states of oscillation are illustrated in Figure 6-2. The 
LCS oscillations were weak, with typical efficiencies of 1.5 percent. 
Although a large 1 2 f 1 1 ratio is associated with high RF efficiency, it 
does not guarantee high efficiency. In one locally optimized repre- 
sentative case noted, we had I /I = 7, but the RF output efficiency 
was less than 0.1 percent, or more than 10 times weaker than with some 
LCS oscillations. Frequency also does not distinguish the two states, 
since the frequency of oscillation has been observed to increase typically 
by only 1.3 percent when a diode spontaneously switches from LCS to 
HCS if Vq or the RF tuning is unchanged. Following application of the 
bias pulse, the LCS oscillations have negligible power but continue to 
grow until some time later (typically 100 ns) as shown in Figure 6-2 (a). 
These oscillations are moderately sensitive to 1^ and are very much less 
sensitive to RF circuit tuning than are HCS oscillations. 

A better-known hallmark of the HCS is that the rapid transition 

from the LCS to the HCS occurs after an interval, t.. We have observed 

t ± to be as low as 40 ns or as high as 300 ns. Lower efficiencies, T], 

and lower Ig/I ratios are generally, but not always, associated with 

longer values of t , whereas high T) and I /I are observed to occur along 

i 2 1 J 
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(a) LOW-CONDUCTIVITY STATE <b) LOW-CONDUCTIVITY STATE SWITCHING 

INTO HIGH-CONDUCTIVITY STATE 


FIGURE 6-2 WAVEFORMS OBSERVED IN PULSED AVALANCHE-DIODE EXPERIMENTS 


with all values of t . For a given R , t. usually first decreases and 

i LL 1 

then increases as V is increased. Otherwise, we have not yet been 

0 

able to identify which properties of the RF circuit or bias circuit 

determine whether t will be large or small. Furthermore, t. is not 

1 1 

constant from pulse to pulse, and the statistical spread in t is 

itself variable. Generally, but not always, the percentage spread is 

low when t is low. Frequently, but not always, t. will increase 
i 1 

substantially, and with no change in 1^, I g , V^, V g , or P when a 
FD-300 diode is chilled with Freon spray. Similarly t may decrease 
when the diode is heated. This suggests that t is related to the time 
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required for the diode to reach a thermal equilibrium consistent with 
the RF oscillations permitted by the circuit. In Section 2 we showed 
theoretically that the period of oscillation is strongly dependent upon 
the magnitude of the diode saturation current. But the saturation 
increases with the increased diode temperature. The experimental obser- 
vation of t being temperature dependent appears to be consistent with 
i 

theory. Additional experiments where the diode is placed in a control- 
led thermal environment are warranted to confirm this point. 

Events occurring during the interval t are of interest. We have 

i 

been able to determine experimentally that the RF oscillations are LCS 
in character and typically, only 1.3 percent lower in frequency than 
that generated in the HCS. These oscillations grow in magnitude until 
the transition to the HCS occurs. Only the first few harmonics of the 
fundamental frequency have been observed, and they were relatively 
small. Whether or not high-frequency "iMPATT" oscillations trigger the 
change of state, as supposed by some, is still uncertain. The con- 
struction of our circuits, however, is such that a high-Q resonance 
at the IMPATT frequency is unlikely. 

It was observed also that the change in states sometimes failed to 
occur on a fraction of the pulses. On occasion this fraction was large. 
When high-efficiency performance from a given diode is accompanied by a 
high percentage of state-switch misses, another RF circuit can some- 
times give the identical performance but without the state-switch misses. 
Investigation of the parameters affecting this instability effect would 
be of value. 


Some of the effects of V and R on oscillator performance have 

0 LL 

been studied. Experimentally, it was observed that the efficiency and 

output power were unaffected by small changes in R . An increase in 

LL 

from 20 to 30 ohms, for example, could be compensated by increases 
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in V , V , and I , and a change in t , but no change in I , V , and the 
oil i £ a 

output power. Larger changes in give different results. An 

increase from R = 20 to 50 ohms or more led to a decrease in the 

LL 

ultimate RF power and efficiency obtainable. Perhaps, during the LCS 

some junction temperature is established that affects the RF efficiency 

after the HCS is reached. With the pulse supply used, the lowest value 

of r =(V - V )/(I -I) obtainable was about 18 ohms; hence we do 

LL 12 2 1 

not yet know whether an optimum value of R less than this exists. 

The effective dc resistivity of the HCS, WaS also • 

noted. For FD-300 diodes, the range 35 < R g < 60 ohms was observed, 

with the best T) obtained for R = 40 ohms. There is a trend for lower 

values of R to be obtained at highar 7], but many exceptions have 
2 

occurred. 


C. Details of Experimental Circuits 
1 . First Diplexer Circuit 

The first circuit used with channel -dropping diplexers had the 
components shown in Figure 6-3. The two diplexers shown schematically 
were designed, constructed, and tested for the frequencies 1500 and 
3000 MHz (fundamental and second harmonic) with a "maximally flat" 
fractional bandwidth of 0.2. A third, analogous, diplexer was prepared 
for the frequency 6000 MHz to be used to control either the fourth 
harmonic in the above case or the second harmonic for a diode having a 
fundamental oscillation at 3000 MHz. 

The diplexers are of a complementary-impedance design, com- 
bining, in series, interdigital bandpass filters and parallel-coupled- 

* 

strip band-stop filters. Adequate selectivity is obtained with two 
reactive elements (n = 2) per filter, and construction is achieved with 
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POINT CONTACT RF 
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SEPARATION AT TWO (or more) FREQUENCIES 





round rods between ground planes and OSM connectors. Computer programs 
were used to select coupled-rod diameters and spacings, and the Hewlett- 
Packard network analyzer was used in tuning up the filters. 

The circuit included a new bias-T designed to present a match 
over a two-decade frequency band extending from well below the lowest 
frequency for which the diode might exhibit negative resistance to well 
above the frequency of the highest harmonic of interest. Thus, in the 
band of each critical frequency the avalanche diode could "see" only 
the load connected to the channel-dropping arm of the respective diplexer, 
while a't all other frequencies it "saw" a "match." This bias-T is com- 
bined with the dummy load for the diplexer and is therefore uncritical 
since it is intended to be dissipative at the outset. The ultra-wide- 
band operation of this bias-T was realized with microstrip conductors 
and physically very small lumped reactive elements. Commercial bias 
networks were not acceptable because of insufficient bandwidth, in- 
adequate peak voltage rating, and excessive capacitance relative to 
the short, fast-rise pulses to be used. 

The system was calibrated initially by replacing the diode by a 
coaxial connector and placing resistive terminations at the 1500- and 
3000-MHz ports. The combined residual VSWR referred to the diode port 
was below 1.3 in the critical ranges 1350-1650 MHz and 2700-3300 MHz. 

In addition, the VSWR was below 1.3 in the low range 100-600 MHz, 
provided that the bias supply looked substantially resistive. The 
VSWR did not exceed 1.9 in the high range 3900-7200 MHz, although the 
1500-MHz diplexer had spurious pass bands 300 MHz wide centered at 
4500 and 7500 MHz. 

The bias-pulse voltage waveform was monitored with conventional 
probes whereas the diode-current pulse envelope was monitored with a 
ferrite-core transformer of novel design. Pulse lengths of 100 to 
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400 ns were used at a PRF of the order of 1 kHz. The bias supply com- 
bined a dc "pre-bias" supply with a low-voltage pulse in a "dc restorer" 


circuit. We now know that this supply was inappropriate for several 
reasons, including inadequate current capacity and lack of resembalance 


to a voltage source with finite R 


LL 


In addition we found that because 


V < V , the diodes exhibited surface breakdowns when a "pre-bias" was 
2 B 

combined with a low-voltage pulse. 


Some Varian experimental TRAPATT diodes were used in the 

circuit configuration shown in Figure 6-3. These diodes typically had 

V « 60 volts and f « 20 CHz. Oscillations were obtained at the 

B IMP ATT 

desired frequency of 1500 MHz although sometimes the fundamental was at 

4500 MHz, which represented a spurious response of the fundamental 

diplexer. The low efficiency and noisiness of the power output suggest 

that the diode was not operating in the HCS. This was further indicated 

by the fact that I /I was less than about 1.2. It is reasonable to 
2 1 

suppose that satisfactory performance would have been obtained with a 
proper bias supply and that useful data could now be obtained with the 
configuration of Figure 6-3. 


2 . Second Diplexer Circuit 

It was decided to work with Fairchild FD-300 diodes, which 

are very inexpensive. From a large batch of these diodes one can select 

diodes having V of any desired value from about 165 V to over 300 V. 

B 

Others have obtained TRAPATT oscillations in the 500-to-1000-MHz range 

13 14 

using these diodes. 2 * * * * * * 9 The diodes are packaged in glass, with thin- 
wire pigtails and an S-shaped contact spring on one side of the silicon 
chip. Thus, these diodes are limited by low thermal dissipation 
capacity both during a pulse and from pulse to pulse. 
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It was also decided to control four harmonic impedances as 
follows. The diode would see adjustable impedances at 750 MHz (n = 1) 
and 2250 MHz (n = 3) through the use of a triplexer and suitable double- 
stub tuners for each of these two channels. The triplexer was built by 
SRI under Contract NAS 12-2165 and is described in the report NASA 
CR-111784. The design is based on a three-reactive-element Chebyshev- 
form fractional bandwidth of 0.286. At 1500 MHz (n = 2) and 3000 MHz 
(n = 4) the diode would "see" an open circuit. This was achieved by 
placing an open-circuited series stub next to the diode. 

The complete arrangement is shown schematically in Figure 6-4. 
Because of the series open-circuited stub it was necessary to design a 
new bias-T that would not seriously load the RF circuit. A schematic 
of the bias circuit is shown in Figure 6-4, including the resistances 
that were added for damping when LC resonances in the bias network were 
observed. The circuit was checked out on the network analyzer and 
found to be satisfactory except for some loss between the diode port and 
the triplexer. Some of this loss was in the stub assembly, which was 
built in microstrip. The rest of the loss,, which was most noticeable 
below 650 MHz, was due to a small coupling to the dissipative elements 
in the bias structure. 

The equipment was operated first with a bias supply now known 
to be unsuitable in a variety of ways, and finally with a supply that 
has been most satisfactory. The pulse current and voltage capacity 
of the bias supply, are adequate for series and parallel combinations 
of FD-300 diodes, and the supply acts as a voltage source with an 
internal resistance of typically 18 ohms. Generally, external series 

resistance was added to optimize R . 

XjLi 

Oscillations were obtained at the expected TRAPATT frequency, 
but they were inefficient and noisy. The oscillations were probably 
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FIGURE 6-4 SECOND EXPERIMENTAL SETUP FOR PULSED OPERATION OF AVALANCHE 
DIODES WITH SEPARATE IMPEDANCE CONTROL (Two fixed, two variable) 

AT FOUR (Harmonic) FREQUENCIES 


100 




in the LCS in that if a 11 jump-up 1 1 in current was observed at all, it 

corresponded to I /I £1.2. Occasionally, we observed a modulation 
2 1 

at about 20 MHz superposed on the pulse waveforms of I , V, and the RF 

envelope of 750 MHz in the region t > t . Similar observations have 

15 1 

been reported by others. Microwave generation was also observed in 
even cheaper Amperex A-13 diodes. These diodes have thicker pigtails 
and no S-spring, but many exhibited excess leakage resistance. Even 
cheaper 1N4454 diodes exhibited low-frequency oscillations. This is a 
new observation not previously reported. 


A possible explanation for the low efficiencies observed was 
suggested from a Smith-chart plot of the impedances looking into the 
diode port of the circuit of Figure 6-4. In the band 643 to 857 MHz 
it was found that there was a small but critical "forbidden zone" in 
the form of a ring on the Smith chart which roughly covered the region 
5 < VSWR < 00 . That is, regardless of the positions of the double-stub 
tuners, the region 0 < Re[z] < 10 ohms was unattainable because of 
circuit losses. The stub tuners were found not to be the principal 
source of loss since measurements indicated that they could give a 
VSWR as high as 25:1 or Re[z] as low as 2 ohms. The losses for the 
fundamental band were partly in the diplexer, but mainly in the series 
stub and bias-T assembly. Recently a new design for the bias circuit 
and circuits at the second and fourth harmonics has been developed which 
should reduce these losses. The circuit of Figure 6-4 with the new 
bias circuit should provide useful knowledge of the TRAPATT mode as 
originally proposed. 


3. Avalanche-Diode Circuit Using Only a Double-Stub Tuner 

Some experiments were made with FD-300 diodes in circuits 
without independent harmonic tuning. If the circuit impedances as seen 


101 




by the diode are readily measurable on a network analyzer, then infor- 
mation on harmonic tuning can be obtained indirectly. The simplicity 
of these circuits permits one to readily gain experience with the high- 
efficiency mode of operation and to measure the diode performance. 

Figure 6-5 shows the diode holder, which used a bypass 
capacitor for the bias input. This method of applying bias has proved 
to be the most useful method for pulse operation where we are not 
limited by thermal dissipation. The bypass capacitance can be as large 
as desired, subject to distortion tolerance on the bias pulse. It is 
possible that this capacitance may affect the speed of the LCS-to-HCS 
transition, Figure 6-2 (b). Our capacitor was formed by clamping a 
0.5-inch X 0.5-inch plate over dielectric tape (typically, 0.001-inch 
mylar) . The corresponding capacitance was roughly 100 to 200 pF, 
depending on the dielectric and the clamping pressure. 

It was useful to measure the impedance looking into the diode 
holder terminals with the diode forward biased sufficiently (typically 
30 mA) to "short out" the junction but not the diode package inductance 
and the holder parasitics. The holder assembly then looked, relative 
to 50-ohm coax, like a low-loss short-circuited stub of length £ , where 
t varied with frequency as given in Figure 6-5. Then if the vector 
impedance of some RF load circuit to the left of the reference plane is 
known, one has only to translate the reference plane through the dis- 
tance 't to refer this impedance to the diode junction. This procedure 
was used to prepare the impedance plots on the following pages. 

The essentials of the rest of the RF circuit are shown in 
Figure 6-6. The procedure was to monitor the pulsed waveforms of I and 
V and the envelope of the fundamental power. The RF signal was detected 
by a calibrated point-contact diode which had a sensitivity close to 
50 mV output across 100 ohms for 0 dBm input. Fundamental and harmonic 
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ELECTRICAL EQUIVALENT OF 
STRUCTURE WITH 1 TO 30 mA 
FORWARD BIAS IN DIODE 





DOUBLE-STUB 

TUNER 

(also dc return) 



FIGURE 6-6 THIRD EXPERIMENTAL SETUP FOR PULSED OPERATION OF AVALANCHE 
DIODES WITH DOUBLE-STUB-TUNER CIRCUIT 

power spectra were also monitored. After many trials, a single 

relatively "superior" diode was chosen, and data were taken for a 

limited number of fixed values of A while optimizing the RF output 

power at the fundamental frequency by varying Y , Y , and V . R was 

12 0 LL 

found not to be critical when kept near its minimum value. Variation 

°f Vg usually yielded only a single maximum whereas all possible 

variation of Y^ and Y^ yielded numerous maxima, which may be referred 

to as local maxima. One of these would be a "grand maximum." V . 

l' 

V 2> 1 g ’ and tuner settings were recorded along with the peak funda- 

mental power output (P^) at the fundamental frequency, f^. Computations 
included R 2 = R T T = (V^ - V^)/ ( 1 ^ ~ and fundamental mode 

efficiency = P^/V^I^. F:i - na lly> the RF circuit was broken at the 


104 





indicated reference plane and the circuit impedances measured with the 
network analyzer at harmonics 1 through 4. 

A selection of representative data (17 cases) is given in 
Figure 6-7 and Table 6-1. In no case was the Re[z] > 20 ohms at any 
of the four harmonics. For n = 1 and n = 2, the Im[z] was almost 
always inductive in contrast to the theoretical results of Section 2 
for the TRAPATT mode. The major exception to the "inductive Im[z]" 
rule is Case G, where tuning is relatively uncritical and Im[Z] is 
substantially capacitive for all four harmonics. In addition, Re[Z] 
tends to be higher than in the other cases, especially at n = 1. Case 
G is included for theoretical interest since it represents the LCS 
and not the HCS despite the same (within 1.5 percent) low frequency of 
oscillation. 

Case D is of some interest because the HCS was clearly obtained, 
but the RF output was even less than can be obtained in the LCS. 
Apparently, wide RF swings of junction current and voltage occurred but 
their relative phases were such that very little power was delivered 
to the load. 

Case L was noted as having the shortest observed t., but this 

feature seems to be un cor related with anything in the impedance plots. 

The same seems true of Case F, in which all three pulse waveforms 

developed a 14-MHz modulation once the HCS was reached. Cases H and 

K exhibited the most power and efficiency although state switching did 

not occur on all bias pulses. The efficiencies obtained compare well 

14 

with those of Chaffin and EerNisse but are somewhat less than was 

13 

typical for Kostishack for single diodes. The impedance-plot 
indications for cases H and K are as follows: 

• At n = 1, relatively high Re[z] among the cases 
presented (4 ohms) 
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Table 6-1 


OPERATING DATA FOR FD-300 AVALANCHE DIODE NO. 54 IN 
DOUBLE-STUB-TUNNER CIRCUIT 
(Case G is LCS, all others are HCS) 


Case 

Fundamental 

Frequency 

Fundamental 
Peak Power 

Fundamental 

Mode 

Efficiency 




Remarks 


f 

1 

(MHz) 

P 

1 

(watts) 

\ 

(percent) 

*2 

(A) 

V 2 

(V) 

R 2 

(0) 


A 

680 

22.4 

5.2 

3.08 

140 

46 


B 

680 

5.0 

1.3 

2.7 

140 

52 


C 

680 

3.2 

0.8 

2.77 

140 

51 


D 

680 

0.25 

0.06 

3.24 

140 

43 

Lower P^ and T\ 
than for LCS 

E 

680 

10.0 

2.6 

2.77 

140 

51 


F 

680 

6.3 pk 


2.3 pk 


61 

Has 14-MHz envelope 
modulation 

G 

670 

0.4 to 1 

0.16 to 0.7 

0.8 

i 

190 

240 

LCS, not HCS (I 2 = 
I 1 , etc.) 

H 

860 

70.0 

14.0 

3.54 

140 

40 

Best efficiency 

I 

865 

1,3 

0.25 

3.7 

140 

38 


J 

860 

10.0 

2.3 

3.16 

140 

44 

Maximum t. 

i 

K 

860 

63.0 

13.0 

3.5 

140 

40 

Best efficiency 

L 

780 

11.0 

2.2 : 

3.5 

140 

40 

Minimum t 

i 

M 

790 

10.0 

2.9 

2.46 

140 

57 


N 

780 

2.0 

0.7 

2.16 

140 

65 


0 

790 

4.5 

1.1 

3.0 

140 

47 


P 

780 

10.0 

2.0 

3.6 

140 

39 


Q 

790 

13.0 

3.3 

2.8 

140 

50 
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• At n = 1, only slightly inductive Im[z] 

• At n = 2, 3, 4, relatively low Re[z] 

• At n = 2, almost no reactance in Z 

• At n = 3, 4, capacitive Im[z] . 

4. Simple Tunable Avalanche-Diode Circuit 

Another experiment was conducted to evaluate further the 

effects of harmonic tuning by indirect means. The RF circuit used is 

13,14 

shown in Figure 6-8 and was related to that used by others. The 

circuit includes a movable short circuit, in a 50-ohm transmission line, 
on one side of the point where the diode is connected to ground via the 
series bypass capacitor. On the other side of the diode, going toward 
the load, a variable, lumped, shunt capacitor was placed at a variable 
location. Thus, three RF tuning variables were available. 

In Kostishack's work the tuning capacitor was about 1/9 
fundamental wavelength away from the diode, and the movable short was 



FIGURE 6-8 ESSENTIAL CONFIGURATION OF “FORTUITOUSLY EFFICIENT" AVALANCHE- 
DIODE CIRCUIT (In Microstrip) 
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located at approximately 1/6 fundamental wavelength from the diode. 

Using the mechanical parameters recorded by Kostishack, and making some 
simplifying assumptions about diode pigtails, we made computations of 
the impedances presented to the diode: 

At 570 MHz, Z i =16.5 + j20 ohms 

At 1140 MHz, Z = 200 + j25 ohms 

At 1710 MHz, Z = 0.05 + j25 ohms 
3 

At 2280 MHz, Z^ = 20 ohms. 

In addition to the very high Re[Z] at n = 2, an inductive Im[z] at n = 

1 and 2 is to be noted. The latter fact is theoretically surprising, 
but it does agree with our experimental observations to date. 

We built a version of this circuit in 50-ohm microstrip, rather 
than coaxial line. Microstrip was chosen primarily for ease in con- 
struction. The sliding short and movable variable capacitor were not 
really suited to microstrip design. The short was imperfect and there 
was inductance in series with the capacitor. These problems were not 
serious, however, for the type of measurements we were interested in 
making. The bypass capacitor was the same as used in the previous 
circuit. With this circuit, we were unable to obtain any stable, 
efficient, HCS oscillations where the dimensions S and T were of the 
order of 0.1 to 0.2 of a fundamental wavelength. Efficient, stable, 

HCS oscillations were observed only when S « T ^ 0.45 of a fundamental 
wavelength. The efficiency compared with the best observed when the 
same diode was used in the circuit of Figure 6-6. The stability was 
improved (i.e. , no state-switch failures on some pulses) and t. was 


jje 

Private communication. 
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generally at a minimum. Typical operating data for FD-300 diode 
No . 54 are : 

f = 850 MHz 

1 

= 50 watts (neglecting radiation) 

P more than 30 dB weaker and P still weaker 
2 3 

t ^ 50 ns 

1 

I = 0.77 A 
1 

I = 3.4 A 

2 

V = 200 V 

1 

V = 140 V 
2 

R = 23 0 
LL 

R = 41 Q 
2 

^ 10.5 percent. 

When another diode (No. 61) having the same V was added in parallel, 
the combination exhibited the best performance with no change in the RF 
tuning adjustments. The fundamental frequency shifted to 825 MHz, t. 
increased to 150 ns, P^ increased by 1 to 2 dB, and T\ increased to 13 
percent. I increased to only 4.2 A, which means each diode had much 
less current. Adding a third comparable diode (No. 30) gave no further 

change other than a slight increase in f . 

1 

The impedances in our circuit were readily measurable on the 
network analyzer. Referred to its junction, the diode "saw" the 
following impedances: 

At n = 1, = 2 + jl3 ohms 

At n = 2, = 7.5 + jlO ohms 
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At n = 3, Z = 12.5 - j25 ohms 

3 

At n = 4, Z = 175 - jlOO ohms. 

4 


5. Avalanche Diodes in Combination 

In addition to what has already been mentioned, some further 

experimentation was done with combinations of FD-300 diodes. These 

data are only preliminary. A series combination of two diodes was tried 

in the circuit of Figure 6-4, but the best performance observed was 

unstable, noisy, and inefficient. Parallel combinations of diodes with 

comparable V 's usually yielded improved stability, and less noise, with 
B 

only a small increase in total I and consequently less current per 

diode. The overall efficiency was generally reduced (except as noted 

previously). It is possible that with a circuit having a high impedance 

at n = 2 the opposite might be true. Kostishack observed increased 

13 

overall efficiency and more noise with diodes in parallel. 

A novel hook-up of two diodes was suggested and investigated 
in part. As shown in Figure 6-9(a), the two diodes are biased in 
series, but to the RF circuit they appear in parallel. If the "DC 
Block” in Figure 6-9(a) is a series open-circuit stub, this circuit 
will achieve the original intentions of the circuit of Figure 6-4 (viz., 
forcing a high impedance to be seen by the diode at certain harmonics 
regardless of tuning further down the line) without requiring any 
troublesome bias chokes or inductors. At present, this is seen as the 
principal application for the diode arrangement of Figure 6-9(a). 

To evaluate this bias scheme, we connected two diodes with 
comparable V^. 1 s as shown in both Figures 6-9(a) and (b) . A commercial 
broad-band "DC Block" was used, having, perhaps, as much as 0.2-dB loss. 
The double-stub tuner was optimized to give maximum output power in the 
fundamental. Tentative data are as follows: 
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LOAD 


LOAD 



(a) BIASING IN SERIES, 
PARALLELING FOR RF 


TA-8327-29 

(b) CONVENTIONAL PARALLELING 


FIGURE 6-9 NOVEL HOOKUP OF TWO AVALANCHE DIODES 


Figure 6-9(a) 

f = 600 MHz 
1 

P = 5 W 

1 

T] = 0.7 percent 

V = 320 V (across 2 diodes) 
2 

I = 2.3 A 
2 


Figure 6-9(b) 


635 MHz 
6 W 

2 percent 
140 V 

2.2 A (2 diodes) 


Remarks: Noisy, with lots of state- Clean, stable HCS ourput 

switch failures 


The data are too incomplete for comment. Future investiga- 
tion should use an improved DC Block or the series open-circuit stub 






mentioned earlier. The diodes should also be selected such that each 

has a similar I operating point for maximum efficiency at a given 
2 

frequency. 


113 



PRECEDING PAGE BLANK NOT FILMED 


7. CONCLUSIONS 

The theoretical understanding of high efficiency oscillations in 

an avalanche diode has been increased by examining three modes of 

oscillation and by developing approximate analytical solutions for 

each mode. A semianalytical solution for TRAPATT oscillations in a 

PIN avalanche diode was obtained by considering seven different phases 

of one cycle and assuming a square-wave current driving the diode. It 

has been shown that a traveling avalanche zone does not form as in the 

P + NN + structure. The resultant voltage waveform is qualitatively 

4 - 4 * 

similar to that obtained for a P NN diode for an assumed square-wave 
current. The highest efficiencies were obtained for the largest drive 
currents and the smallest saturation currents. In many instances the 
square-wave drive current does not correspond to the diode being 
embedded in a passive circuit. 

Based on the analytical solution, an approximate computer model of 
TRAPATT oscillations in a PIN avalanche diode was developed. This 
model has a computer running time at least two orders of magnitude less 
than more exact computer models. The model was used to calculate the 
diode voltage waveform and operating characteristics for several 
different assumed diode current waveforms. A square-wave current was 
approximated by using the first few terms of its Fourier series. It 
was found that both the period of oscillation and the dc-to-RF con- 
version efficiency decreased from those of a corresponding square wave. 
The higher harmonics are necessary for high efficiency because they 
allow the current to rise rapidly and drive the diode far into avalanche 
before appreciable charge is generated. Subsequently a dense plasma 
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is formed , forcing the diode voltage to be very low while the current 
remains high. Equally important for high efficiency is that the 
current be zero or nearly zero over the second half of a period so 
that the diode voltage remains high over this half of the cycle. 

There are numerous voltage and current waveforms consistent with 
the device physics of an avalanche diode that corresponds to high- 
efficiency operation at the fundamental frequency. Not all such sets 
of waveforms, however, correspond to the device being embedded in a 
passive microwave circuit with only a dc source. The results of 
Section 2 show, qualitatively at least, that high-efficiency oscilla- 
tions are possible with currents containing only one or two harmonics. 

A second mode of high efficiency oscillations which is very similar 
to the ordinary TRAPATT mode was described and a semianaly tical solu- 
tion developed for the voltage for a rectangular current waveform. 
Portions of a cycle of oscillation were shown to be identical with those 
of the ordinary TRAPATT mode, but during avalanche the external current 
was assumed to be small. No calculations have been made for this mode 
although it is believed that it may correspond to experimental obser- 
vations of moderately high and high efficiency oscillations. 

Another significant result obtained in this report is the develop- 
ment of a model of high efficiency oscillations in a PIN diode for 
extremely large current densities which causes avalanche to occur over 
a major portion of a half cycle. This mode has not been discussed in 
the literature before. We showed that the output power of the device 
increases proportionally with the magnitude of the current with only a 
relatively small decrease in dc to rf conversion efficiency. In 
practice operation in this mode would require circuits with lower 
impedances and better heat dissipation. 
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An approximate solution for the voltage waveform in a PIN diode 
was also derived using an idealized step function for the ionization 
coefficient. The utility of this approximate model is that it can be 
used to indicate possible patterns of avalanche regions in the space- 
time plane. Analytical solutions can then be written for each region 
and the solutions matched at the boundaries between the avalanching 
and nonavalanching regions. Estimates of the current and voltage 
variation across a diode were deduced without computation for one 
permissible pattern. It was shown that within limitations the current 
waveform could be selected for efficient device operation. 

Experimental observations of high efficiency oscillations were 

made using inexpensive diodes. The highest efficiencies and output 

powers were comparable to those reported by others for these diodes 

13 

except for the anomalous results reported by Kostishak. Simpler 
circuits than those originally proposed for independent harmonic tuning 
were used. The efficiencies obtained with the channel dropping filters 
for harmonic tuning were typically less than a few percent. These 
results, however, are inclusive because the circuit losses were too 
large and the bias supply inadequate for these early measurements. A 
disadvantage of using filters to separate the harmonics for independent 
tuning is that the tuning elements are a long distance away from the 
diode. The tuning becomes more critical and the losses are increased. 
Ideally, impedance matching for each harmonic should be accomplished 
next to the diode for optimum efficiency. 

The measured circuit impedances up to at least four harmonics of 
those circuits which gave highest efficiencies did not correspond to 
the theoretical results for TRAPATT oscillations obtained in Section 2 
of the report. This can be explained by the fact that many different 
voltage and current waveforms are possible; all of these lead to high 
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efficiency operation. In our analyses we assumed current waveforms 
composed of a few harmonics and calculated the corresponding voltage 
waveform which always had several higher harmonics. An equally valid 
method would be to assume a voltage waveform composed of one or a few 
harmonics and calculate the corresponding current waveform which may 
be rich in harmonic content. The latter corresponds to circuits with 
very low impedances at the higher harmonics. This is more representa- 
tive of the experimental observations. On the other hand, analysis of 
13 

Kostishack's circuit indicates that the impedances at the even 
harmonics are large, which corresponds more closely with the theoretical 
results of Section 2. 
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8. RECOMMENDATIONS 


In as much as the experiments with tuning at multiple frequencies 
were inconclusive, it is recommended that these be repeated with the 
improved bias supply now available. It is also recommended that simpler 
circuits be designed which take into account the desired impedances 
at some of the higher harmonics but do not provide flexibility in tuning. 
The impedance matching could be designed to be close to the diode and 
the loss could be kept small in these circuits. By building a limited 
number of these circuits, one could verify the theoretical results 
obtained to date. 

In the course of the experimental work several observations were 
made which led to some interesting theoretical and experimental questions. 
On any one bias pulse the time that the diode was in a low conduction 
state before switching into high efficiency oscillations appeared to be 
temperature dependent. For pulse applications this delay time would be 
important. Controlled experiments should help to identify the critical 
parameters affecting this delay time. 

It was also noticed that from bias pulse to bias pulse the diode 
did not always begin to oscillate at the same time. In fact, on some 
pulses the diode did not break into oscillation at all. These effects 
may also be thermal in nature. A significant problem in the high 
efficient modes of operation is the noisy nature of the oscillations. 

It is recommended that a significant theoretical and experimental effort 
be directed toward understanding the noise properties of high efficiency 
oscillations in avalanche diodes. 
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Although the theoretical results obtained in this report are 
significant, they are not complete. First, more extensive calculations 
should be made of the second high efficiency mode and the mode with 
multiple avalanche discussed in this report. These calculations could 
be compared with the ordinary TRAPATT mode and guidelines developed for 
optimum operation and corresponding circuit design. The simple but 
economical computer program for TRAPATT oscillations in a PIN diode 
should be extended to include other doping distributions and nonequal 
ionization coefficients and carrier velocities. A similar economical 
computer program should be developed to calculate the diode current when 
the voltage waveform is specified. 

Accurate and detailed computation of the operating parameters of 
an avalanche diode oscillating in a high efficiency mode require a 
complex computer program like UNSAT developed for this contract. However, 
the running time is extensive when the charge carrier and field gradients 
are large as they are in the TRAPATT mode. However these gradients 
usually exist only over a small portion of the diode depletion region. 

It may be possible by more elaborate programming to incorporate a 
variable step size in the spatial coordinate. Where the gradients are 
large the step size would be small and where they are small the step 
size would be large. This technique would permit accurate calculations 
more economically. 
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Appendix A 


NORMALIZED VARIABLES 


Below is the normalization used and the values each variable should 
be multiplied by to obtain unnormalized values. 


Unit 

Expression 

Value 

Length 

w 

-3 

10 cm 

Time 

w/v 

-10 

10 second 

Charge Density 

s 

N 

o 

-15 3 

10 cm 

Current 

qv N 
so 

3 2 

1 . 6 X 10 A/ cm 

Field 

wqN 0 / e 

5 

1.6 x 10 volts/cm 

Unit Charge 

q 

-19 

1.6 X 10 coulomb 

Voltage 

2 

w qN Q /e 

160 volts 

Impedance 

2/ 

w /e v 
s 

2 

0.1 ohms -cm 

Power 

(qwN 0 ) 2 v s /e 

5 , 2 

2.56 X 10 watts/cm 
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Appendix B 

AN ECONOMICAL COMPUTER PROGRAM FOR TRAP ATT OSCILLATIONS 

In this appendix we describe an approximate time-domain model for 
TRAPATT oscillations in a PIN avalanche diode. During these portions of 
a cycle when the E-field is effectively uniform, Eqs. (2-8) and (2-10) 
of Section 2 express the dynamics of the total internal charge and the 
terminal voltage of the diode. When a trapped plasma has formed and is 
being removed, the electric field increases from E p in the plasma to a 
larger value at the edges of the depletion layer because of the residual 
charge. Let x be the distance from the right-hand edge of the depletion 
layer to the corresponding edge of the trapped plasma. It is also the 
distance from the left-hand edge to the plasma. Then, during the field 
depression phase of the TRAPATT cycle, x and the diode terminal voltage 
are given by the following two equations. 


dx 

dt V 


(B-l) 


2 

V = n x + E 
r P 


(B-2) 


To a first approximation the residual charge, n^, is equal to J/v^ ^ 
the time variations are not too rapid. This has been verified by more 
accurate computer simulations using program UNSAT. Thus, Eq. (B-2) 
becomes 

V(t) = X 2 + E (t) . (B-3) 

V p 

s 
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During the plasma extraction phase, the edges of the plasma move at a 

plasma velocity, v ; thus Eqs. (B-l) and (B-3) become 

P 


dx 

dt 



(B-4) 


✓ x J(t) 2 , v 

V(t) = x + E (t) , (B-5) 

v p 9 

s 

E 

where E (t) is obtained by inverting the expression for v = 

P p E+0.1 ’ 

Equations (B-l) through (B-5) are valid when Q > J/v^. When the charge 
density is so small that the particle current (at a saturated velocity) 
is less than the external current, the E-field is nearly uniform across 
the depletion layer, and Eqs. (2-17) and (2-19) are applicable. 

A computer program (Called modified SMPAL) was written which 
integrates Eqs. (2-17) and (2-19) numerically and uses Eqs. (B-l) and 
(B-3) or (B-4) and (B-5) as appropriate during a cycle. Note that in 
this formulation the more exact expression for the ionization coefficient 
can be used as well as an analytic expression for the velocity as a 
function of the electric field. 

The validity of modified SMPAL was verified by comparing the results 
with a more exact computer simulation, UNSAT. Figure B-l shows the 
charge generation, G, as computed by both UNSAT and modified SMPAL 
assuming that the diode current was constant. The generation predicted 
b y ^be modified SMPAL is too large. This is a consequence of assuming 
that the charges drifting out of the depletion layer can be modeled by 
an exponential decay. In the more exact solution, UNSAT, when the diode 
is biased below avalanche, there is some very small but finite charge 
generated which increases the charge density over that corresponding to 
the saturation current. On the other hand, in modified SMPAL the 
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NORMALIZED SATURATION CURRENT — J s 


TA-8327-13 

FIGURE B-1 TOTAL CHARGE GENERATED IN A PIN DIODE AS 
COMPUTED BY THE COMPUTER PROGRAMS UNSAT 
AND MODIFIED SMPAL FOR A CONSTANT CURRENT 






charge decays when a is less than unity. Therefore, the charge in 
modified SMPAL will initially decrease from the saturation level when 
we start with the E-field below avalanche. With a smaller charge density 
(equivalently saturation current) the diode eventually charges up to a 
higher voltage, increasing the total charge generated. We can compen- 
sate for the differences in generated charge in the two programs by using 

a larger value of J in modified SMPAL than the corresponding value in 

s 

program UNSAT. By doing this we are able to compare the voltage wave- 
forms computed by both programs during the initial charge-up and genera- 
tion phases of a TRAP ATT cycle. Figure B-2 shows the voltage waveforms 
during the generation and field depression phases of a TRAPATT cycle as 
computed by modified SMPAL and UNSAT for two different current levels. 

The saturation currents used in each program were obtained by choosing 
the same value of generation for each and using the corresponding values 




TRANSIT TIMES AT SATURATED VELOCITY TRANSIT TIMES AT SATURATED VELOCITY 


TA-8327-12 


FIGURE B-2 COMPARISON OF THE TERMINAL VOLTAGE OF A PIN DIODE DURING THE 
INITIAL CHARGE-UP AND GENERATION PHASES OF A TRAPATT CYCLE 
AS COMPUTED BY PROGRAMS UNSAT AND MODIFIED SMPAL 
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of J given by the respective curves in Figure B-l . The agreement 
s 

between the results from the two programs is excellent. The waveforms 
were not compared over an entire TRAPATT cycle because the exact 
solution was too expensive to run. 
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